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Abstract  (cont'd) 

tunnel  traversing  mechanism.  The  associated  vorticity  field  was  inferred  from 
these  measurements  through  a spatial  contour  integration  procedure.  < 

The  existence  of  several  identifiable  "flow  regions"  exhibiting  similar 
characteristics  was  established  and  verified  with  both  hot-wire  and  flow 
visualization  data.  Among  the  regions  examined  extensively  are  the  three- 
dimensional  wing  boundary  layer,  the  inboard  potential  flow  region  and  tne 
region  of  vortex  shedding  and  roll-up  adjacent  to  the  wing  tip. 

Measurements  of  the  three-dimensional  boundary  layer  were  employed  to  deter- 
mine the  accuracy  of  the  measurement  technique  close  to  the  solid  surface  boun- 
dary and  to  insure  that  the  boundary  layers  over  the  wing  remained  thin  and 

~ i^The  inboard  potential  flow  region  was  observed  to  exhibit  the  simplifying 
cliaracteristic  of  a spanwise  cross-flow  velocity  component  which  is  independent 
of  the  surface  normal  coordinate.  As  a result,  a possible  ambiguity  in  the 
determination  of  the  stream-wise  vorticity  component  was  removed.  Graphs  of 
the  near^surface  (bound)  vorticity  distribution  are  presented  along  with  a 
map  of  rfPrandtl  bound  vortex  filaments. "‘"'’The  associated  bound  circulation 
function  could  be  expressed  approximately  as  the  product  of  two  separable 
functions  of  the  planform  spatial  coordinates.  Measurements  at  the  wing 
trailing  edge  are  employed  to  show  that  the  properly  non-dimensionalized  char- 
acteristic vorticity  distribution  in  this  region  is  independent  of  angle  of 

attack  .<1— . . 

Shed  vorticity  emanating  from  the  lower  aerodynamic  surface  near  the  wing 
tip  was  observed  to  roll-up  adjacent  to  the  tip  and  roll  over  onto  the  upper 
wing  surface  at  a predictable  chordwise  location.  The  presence  and  behavior 
f this  structure  which  possesses  an  identifiable  "core"  is  seen  as  an  explana- 
)n  to  the  modified  pressure  distribution  (lift  increment)  normally  observed 

in  this  region.  * . 1 

A comparison  of  the  flowfield  in  the  near  wake  region  with  this  data  reveals 
that  roll-up  occurs  much  more  violently  near  the  solid  surface  than^ in  the 
boundary  free  zone  downstream  of  the  wing.  The  normalized  description  of  the 
roll-up  process  in  the  near  wake  region  is  observed  to  be  nearly  independent 
of  both  aspect  ratio  and  angle-of-attack. 
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NOMENCLATURE 


A - Area,  usually  referring  to  that  enclosed  by  a 

contour 


. , A(x)  Function  expressing  chordwise  variation  of 

circulation 

AR  - Aspect  ratio 

B - A constant  in  the  hot-wire  response  equation 

B(y)  - Function  expressing  spanwise  variation  of 
circulation 

CD  - Wing  drag  coefficient 

CL  - Wing  lift  coefficient 

D - Drag,  usually  out  of  the  wing 

E - Output  voltage  of  hot-wire  bridge 

E - Output  voltage  of  hot-wire  bridge  at  zero 

velocity 

E - Average  output  voltage  of  hot-wire  bridge  (Note: 

overbar  denotes  "average";  usually  with  respect 
to  time.  If  the  averaged  quantity  appears  with 
an  adjacent  variable  in  parentheses,  the  average 
is  taken  with  respect  to  that  variable,  e.g., 

I’  (x)  is  an  average  with  respect  to  the  x coor- 
dinate.) 

I . - Contribution  to  "A"  type  contour  integrals  by 

s flow  field  in  the  vicinity  of  the  leading  edge 

stagnation  line 

I . - Contribution  to  "A"  type  contour  integrals  by 

j 1 flow  field  over  that  portion  of  the  contour 

intercepting  the  boundary  (or  viscous)  layer 

i K - A constant,  sometimes  subscripted  if  necessary 

L - Lift,  usually  of  the  wing 

M - Slope  of  the  lift  curve,  CL  = Ma 

N - A dimensionless  parameter  relating  the  ratio 

of  crossflow  velocity  to  mainflow  velocity 


V 


P(x.)  - Integral  of  the  pressure  differential  over  the 

local  span, 

£DAp(x,y)dy 

2 

Q - Dynamic  pressure  of  the  free  stream,  1/2  pUQ 

S - Surface  area,  usually  of  the  wing 

U - General  symbol  for  velocity,  subscripted  to 

denote  specific  spatial  components  (Ux,  U , 

U ) . Also  used  as  velocity  component  y 
afong  hot-wire  probe  direction 

UQ  - Free  stream  tunnel  velocity 

Uii  - Component  of  velocity  in  mainflow  direction 
I'  parallel  (locally)  to  the  wing  surface 

U.  - Azimuthal  component  of  velocity  in  cylindrical 
coordinate  system 

V - Velocity  component  perpendicular  to  hot-wire 

prove  direction 

W - Velocity  component  perpendicular  to  plane  of 

X-probe 

b - Wing  semispan,  geometric 

b*  - Effective  wing  semispan,  with  respect  to  the 
spacing  of  trailing  vortices 

c - Wing  chord  length,  geometric;  or  symbol  for 

"contour" 

c^  - Wing  section  lift  coefficient 

e'  - /Fluctuating,  component  of  hot-wire  bridge  output 

voltage,  E = E + e',  prime  denotes  a fluctuating 
quantity 

h - A thickness,  or  thickness  parameter — used  as 

slab  thickness  or  that  of  body  plus  displacement 
regions 

1 - path  distance  along  a contour,  c. 

1*  '•  Distance  from  wing  trailing  edge  to  a point 

downstream  where  vortex  has  rolled  up 

o - Symbol  meaning  "of  the  order" 
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p - Pressure 

- Pressure  on  lower  (aerodynamically)  wing  surface 

Pu  - Pressure  on  upper  (aerodynamically)  wing  surface 

r - Radial  co-ordinate  in  cylindrical  coordinate 

system;  also  radius  of  curvature 

t - Wing  thickness,  usually  a function  of  x,  t(x) 

u'  - Fluctuating  component  of  .velocity  in  x-direction, 

along  probe  direction  of  hot-wire  probe 

v'  - Fluctuating  component  of  velocity  perpendicular 

to  hot-wire  probe  direction 

x - Flow  direction  parallel  to  free  stream 

y - Flow  direction  perpendicular  to  free  stream 

extending  spanwise  along  the  wing 

z - Flow  direction  perpendicular  to  free  stream 

extending  normal  to  x-v  plane;  also  the  axial 
component  in  cylindrical  coordinates 

a - Wing  angle  of  attack;  also  hot-wire  crossflow 

angle 

T - Circulation 

To  - Circulation  of  a contour  enveloping  the  wing 

plus  boundary  layer  at  the  centerline  spanwite 
location 

T - Circulation  around  a contour  oriented  in  the  xz- 

plane  which  surrounds  all  or  part  of  the  wing 

rB  - Circulation  around  a contour  oriented  in  the  yz- 

plane  which  surrounds  all  or  part  of  the  wing 

A - "The  change  in"--when  used  as  a prefix  to  another 

symbol  it  implies  the  difference  in  two  values 
of  that  quantity.  In  this  paper  used  exclu- 
sively to  reflect  the  difference  between  values 
on  the  upper  and  lower  surfaces  for  the  same 
chordwise  and  spanwise  locations 

6 - Boundary  layer  thickness,  sometimes  subscripted 

to  denote  the  reference  surface 

t,  - Vorticity 
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X - Vorticity  vector,  fk=£'  i + £j  + £k  ! 

x y z I 

A - Wing  tip  sweep  angle  i 

j 

p - Density  of  fluid,  usually  air 

i ! 

a - Sensitivity  coefficient  of  hot-wire  anemometer  j 

4>  - X -me ter  wire  inclination  angle  I 

1 * 1 

0 - Complement  of  X-meter  wire  inclination  angle, 

<J»  = ir/2  - 0;  also  used  as  azimuthal  coordinate 

when  referring  to  a cylindrical  coordinate  - j 

system  j 


CHAPTER  I 


INTRODUCTION 

This  investigation  is  concerned  with  the 
formation  and  early  development  of  a trailing  vortex 
resulting  from  lift  generation  on  a finite  wing.  The 
process  is  examined  experimentally  employing  measure- 
ments of  the  velocity  and  vorticity  fields  adjacent  to 
and  immediately  downstream  of  a rectangular  wing 
immersed  in  a uniform  stream.  The  experiments  were  con- 
ducted in  a low  speed  wind  tunnel  and  involve  nominally 
unseparated  flow' conditions  over  the  lifting  face. 

The  significance  of  the  problem  can  clearly 
understood  by  considering  the  entire  flow  field  which 
results  from  the  passage  of  fluid  over  a lifting 
surface  (Figure  1) . The  disturbances  resulting  from 
the  wing's  presence  lead  to  the  appearance  of  a wake 
which  persists  for  a long  time  period  (far  downstream 
of  the  wing)  and  which  evolves  through  several  identi- 
fiable stages. 

The  flow  near  the  solid  surface  is  classically 
termed  the  "boundary  layer" — that  region  which  is  the 
source  of  all  vorticity  generated  in  the  flow  field. 

The  behavior  of  the  flow  in  this  region  together  with 
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the  neighboring  "external"  flow  is  directly  related  to 
the  surface  load  distribution  and  is  the  primary  arena 
of  interest  in  this  study.  This  region  is  more  fully 
discussed  below. 

The  sector  immediately  downstream  of  the  wing, 
referred  to  as  the  "near  wake"  or  rollup  region,  is 
usually  strongly  three  dimensional.  Vortex  filaments 
shed  into  this  area  are  free  to  influence  each  other 
which  results  in  the  ultimate  development  of  two  dis- 
crete vortices.  This  region  is  one  of  considerable 
importance  in  aircraft  applications  since  a good  portion 
of  an  aircraft  structure  (tail  surface  and  controls)  is 
usually  immersed  in  it.  The  behavior  of  this  flow  can  also 
influence  the  behavior  of  other  aircraft  flying  in 
this  region,  e.g.,  formation  flying,  aerial  refueling. 

The  mechanics  of  rollup  is  the  subject  of  a number  of 

recent  investigations  most  of  which  are  numerical  in 
1-5 

nature.  It  is  examined  to  a limited  extent  in  the 

present  work. 

Further  downstream  the  wake  is  observed  to  be 
composed  of  two  separate  vortical  structures  having 
opposite  directions  of  rotation.  Each  of  these  highly 
stable  forms  is  characterized  by  a cylindrically  sym- 
metric shape  and  composed  of  a viscous  core  and  an  outer 
flow  which  behaves  like  that  induced  by  a potential  line 
vortex.  Frequently,  an  axial  deficit  or  excess  is 
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detected  in  the  core  region  which  prohibits  its  being 
analyzed  as  a two  dimensional  field.  The  behavior  of 
this  part  of  the  flow  field  has  drawn  increased  attention 
recently  due  to  the  hazard  it  poses  to  encountering 
aircraft.  Its  persistence  in  high  density  airport 
traffic  areas  has  been  the  cause  for  much  concern 
especially  when  generated  by  heavy  "jumbo- jet"  type 
aircraft. 

A host  of  investigators  have  studied  the 

6-9 

motion  of  these  vortices,  including  their  stabil- 

ity,'1'0-'1'2 their  interaction  with  various  environmental 

12  13  14  15 

conditions,  ' and  their  decay.  ' Prediction 

methods  for  their  interaction  with  encountering  air- 

1 6 

craft  have  also  been  developed.  They  have  been 

17  18  19-21 

studied  theoretically,  ' numerically  and  exper- 

22-26 

imentally  in  both  laboratory  experiments  and  in 

27  28 

flight  tests.  ' In  short,  many  aspects  of  the  flow 
in  the  wake  region  are  reasonably  well  documented. 

This  is  not  the  case,  however,  with  the  flow  in 
the  immediate  vicinity  of  the  wing.  Though  the  problem 
of  lift  generation  was  one  of  the  first  addressed  by 
fluid  mechanists  at  the  beginning  of  this  century,  the 

t 

primary  objective  was  and  continues  to  be  the  determin- 
ation of  airloads  on  the  wing  motivated  primarily  by 
structural  and  design  considerations.  The  qualitative 
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nature  of  the  flow  was  sketched  by  Lanchester,  and 
. 30 

visualized  by  Hoerner,  but  little  has  since  been  done 
to  ascertain  its  detailed  character  in  the  immediate 
vicinity  of  the  wing. 

Early  experimental  efforts  concentrated  largely 

on  the  measurements  of  integral  effects  (e.g.,  lift 

coefficient)  and  the  surface  pressure  distribution. 

Direct  information  regarding  the  flow  field  was  obtained 

primarily  through  simple  flow  visualization  techniques. 

Analytical  studies  were  confined  to  invisid 

models  since  the  pressure  distribution  was  found  to  vary 

only  slightly  through  the  boundary  layer  over  a large 

portion  of  the  wing  surface.  In  fact,  the  qualitative 

aspects  of  the  flow  field  were  usually  prescribed 

a priori,  reference  Prandtl ' s lifting  line  theory. 

31  32 

Recent  numerical  efforts  ' which  have  attempted  to 
predict  the  structure  of  the  vorticity  field  are  ham- 
pered by  the  necessity  for  user  supplied  "initial 
conditions"  for  the  flow.  An  iteration  scheme  is  usually 
employed  in  which  the  flow  field  is  allowed  to  deform 
slightly  which  in  turn  leads  to  a new  set  of  interaction 
conditions.  The  procedure  is  repeated  until  a "stable" 
situation  is  attained  depending  on  the  operating  limits 
of  the  computer  employed. 

Another  shortcoming  of  most  finite  wing  analyses 
is  the  calculated  behavior  of  the  flow  near  the  wing 


tip.  A mathematical  singularity  in  the  vorticity 

distribution  corresponding  to  an  infinite  velocity  is 

usually  found  there  which  is  physically  impossible. 

Difficulties  encountered  in  further  calculations  are 

usually  resolved  via  mathematical  methods  which  remove  or 

33 

ignore  the  singularity. 

The  actual  physical  situation  is  related  to 
the  analytical  model  through  the  concept  of  bound  vor- 
ticity. In  potential  theory  this  entity  is  usually 
ascribed  to  the  "body"  in  contrast  with  free  vorticity 
which  remains  with  the  fluid  particles  (or  elements) . 

In  reality,  the  bound  vorticity  is  actually  that  which 
is  continually  generated  by  shear  in  the  wing  boundary 
layer  and  which  appears  "attached"  to  the  wing  in  the 
appropriate  Eulerian  reference  frame.  In  two-dimensional 
flow  or  at  the  symmetry  plane  of  a finite  wing, the 
vorticity  vector  is  perpendicular  to  the  free  stream 
direction. 

Of  equal  significance  is  the  apparent  redirec- 
tion of  vorticity  toward  the  free  stream  direction, 
i.e.,  the  horseshoe  vortex  concept.  The  actual  re- 
orientation of  vorticity  is  caused  by  the  spanwise 
variation  of  the  pressure  gradient  over  the  surface. 

If  one  considers  the  effects  on  a fluid  element  which 
lies  near  the  surface  as  shown  in  Figure  2,  it  is 
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observed  to  deform  and  rotate  due  to  the  influence  of 
the  crossflow  velocity  gradient  as  well  as  that  in  the 
mainflow  direction.  As  it  moves  within  the  boundary 
layer  it  continues  to  be  effected  by  this  rotational 
tendency  depending,  of  course,  on  the  strength  and 
character  of  the  crossflow  shear.  Eventually  the 
element  enters  the  free  field  at  some  point  of  boundary 
layer  separation  (trailing  edge,  wing  tip  or  other 
location  where  the  layer  separates)  where  it  becomes 
part  of  a free  vortex  filament  in  the  wake.  The  above 
description  is  not  strictly  accurate  since  the  hypo- 
thetical element  is,  of  course,  distorted  and  diffused 
as  it  moves  over, even  a small  distance,  but  the  discus- 
sion serves  a pedagogical  purpose  in  describing  the 
origins  of  mainstream  directed  vorticity. 

In  addition,  this  simplistic  view  is  compli- 
cated by  the  presence  of  free  vorticity  near  the  sur- 
face which  is  not  accounted'  for  in  the  "sheet"  models 
employed  to  describe  most  wings.  This  vorticity,  like 
that  shed  at  the  wing  tip,  can  alter  the  predicted 
pressure  gradient  which  can,  in  turn,  reinfluence  the 
flow  geometry.  Hence  a simple  cause-effect  relation- 
ship can  never  be  established. 

These  difficulties  are  discussed  in  the  ensuing 
pages  in  an  attempt  to  provide  insight  into  the  physical 
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mechanisms  operating  on  this  type  of  flow  field. 
Hopefully, the  experimental  evidence  uncovered  in  this 
investigation  will  serve  to  aid  in  describing  flows 
over  other  wing  configurations. 


CHAPTER  II 


Most  experiments  were  conducted  in  the  2 foot 
x 2 foot  cross-section  low  speed  wind  tunnel  at  the 
University  of  Colorado  (see  Figure  3).  This  open  return, 
suction  type  tunnel  is  driven  by  a single  stage  axial 
flow  impeller  powered  by  a variable  speed,  direct  cur- 
rent motor.  The  useful  velocity  range  of  this  facility 
extended  from  5 to  110  feet  per  second. 

Air  drawn  into  the  settling  chamber  was  initially 
filtered  through  a furnace  filter  type  dust  trap  at  the 
inlet.  After  passing  through  a hexagonal-honeycomb 
structure  and  several  screens  it  was  accelerated  into 
the  test  section  through  a 16:1  area  ratio  contraction. 

The  resulting  free  stream  turbulence  level 

"IT  1/2 

( (?u.  ) ' /U  ) was  measured  at  0.06%  for  U =60  feet  per 
i i o o 

second.  The  air  was  then  exhausted  back  into  the  room 
through  a conical  diffuser. 

The  14  foot  long  test  section  is  composed  of 
two  separate  units . The  first  unit  (nearest  the  inlet) 
is  4 feet  in  length  and  constructed  of  transite 
insulating  material  on  all  sides.  A small  quartz  window 
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is  located  on  one  wall  as  a provision  for  optical  heat- 
ing of  immersed  bodies.  The  remaining  10  foot  section 
is  constructed  of  plywood  with  a formica  finish  on  three 
sides.  The  fourth  wall  is  a solid,  1/2  inch  thick 
plexiglas  sheet.  This  latter  section  was  the  one 
utilized  during  the  course  of  this  experiment  primarily 
due  to  the  need  for  optical  viewing,  which  will  be  dis- 
cussed below.  Variation  of  the  static  pressure  gradient 
was  less  than  1%  over  the  entire  test  section  length. 
Cross-sectional  variation  of  mean  velocity  was  less  than 
0.4%  outside  the  boundary  layer  in  the  empty  tunnel. 

The  tunnel  environment  was  a large  air-condi- 
tioned room  nominally  maintained  at  70°F.  Thermal 
equilibrium  with  the  tunnel  operating  (60  ft.  per  sec.) 
was  reached  in  approximately  two  hours.  All  data  were 
obtained  with  the  system  in  total  thermal  equilibrium. 

A precision  traversing  mechanism  was  constructed 
to  provide  accurate  positioning  of  the  measuring  probes 
(see  Figure  4 for  a photograph  of  this  mechanism)..  The 
center  body  could  be  moved  either  vertically  or  horizon- 
tally with  a drive  provided  by  miniature  d.c.  motors. 
Movement  was  accomplished  through  a screw  drive  for 
vertical  movement  and  a rack-and-pinion  system  for  hori- 
zontal movement.  Geared  potentiometers  were  utilized  to 
monitor  distances  in  both  these  directions.  The 
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resultant  system  spatial  resolution  was  ±0.005  inches 
in  horizontal  movement  and  ±0.002  inches  in  vertical 
movement.  Axial  movement  of  the  entire  traversing 
system  was  accomplished  manually  by  adjustment  of  a 
connecting  cable  which  runs  the  length  of  the  usable 
test  section. 

A Gaertner  optical  cathotometer  was  used  to 
establish  reference  distances  between  the  model  and 
the  probes  in  various  configurations.  Utilizing  this 
device,  probes  could  be  positioned  to  within  .005  inches 
of  the  model  surface  without  establishing  physical 
contact.  In  addition,  the  instrument  was  used  to 
monitor  the  calibration  of  the  vertical  positioning 
system  of  the  traversing  mechanism. 

Flow  visualization  photographs  were  obtained  in 
the  2 foot  x 3 foot  cross-section  subsonic  wind  tunnel 
at  the  USAF  Academy.  This  closed  circuit  facility, 
powered  by  a motor-generator  driven  single  stage  fan  has 
a useful  test  speed  range  of  50  to  400  feet  per  second. 

A Sage  Action,  Inc.,  soap  bubble  generator  was  located 
in  the  stilling  chamber  just  upstream  of  the  test 
section.  Helium  filled,  neutrally  buoyant  soap  bubbles 
were  created  and  injected  into  the  free  stream  at 
approximately  the  free  stream  velocity.  These  bubbles 
were  then  convected  into  the  test  section  * id  utilized 


i 


11 

as  the  medium  for  observing  flow  field  behavior.  The 
plexiglass  test  section  side  walls  provided  a window  for 
photographic  purposes. 

II. 2.  The  Test  Wing 

The  selection  of  a suitable  lifting  shape  for 
the  experiment  was  predicated  on  both  aerodynamic  con- 
siderations and  facility  limitations.  A simple  geometry 
having  a thin  cross-section  seemed  desirable  in  order  to 
simplify  the  theoretical  aspects  of  the  problem.  An 
airfoil  shape  which  reasonably  satisfied  the  conditions 
of  thin  airfoil  theory  while  yielding  minimal  wake  and 
boundary  layer  displacement  characteristics  was  required 
to  physically  emulate  the  simple,  potential  flow  theo- 
retical models.  Since  the  detailed  flow  field  about  the 
wing  was  to  be  investigated,  considerations  of  the  posi- 
tion resolution  capabilities  of  the  apparatus  and  the 
minimum  effective  probe  size  demanded  a reasonable  wing 
chord.  The  relatively  small  tunnel  cross-section, 
however,  imposed  a severe  restriction  on  the  span  para- 
meter thus  limiting  the  maximum  effective  aspect  ratio. 
(Initial  measurements  behind  a wall-mounted  rectangular 
wing,  NACA  2412  cross-section,  with  half-aspect  ratio 
3 indicated  that  almost  all  vorticity  production  occurred 
within  a 1 chord  length  distance  of  the  tip  in  the  span- 
wise  direction.) 


12 


On  these  bases  a rectangular,  untwisted  wing 
having  an  NACA  64009  laminar  flow  airfoil  cross- 
section  was  constructed  from  solid  wood  with  a sharp, 
aluminum  trailing  edge.  The  surface  was  finely  polished 
in  an  attempt  to  preserve  the  laminar  flow  characteris- 
tics. Maximum  surface  ripple  was  less  than  0.02  inches 
over  the  entire  span  of  the  model . The  surface  rough- 
ness was  minimal'  except  for  two  small  irregularities  on 
the  inboard  portion  of  the  leading  edge.  To  provide  an 
accurate  geometric  reference  for  the  measurements,  dark 
black  spanwise  lines  were  drawn  over  the  white  finish  on 
both  upper  and  lower  surfaces  and  extended  around  the 
exposed  tip  side  at  5%  chord  increments  from  the  leading 
to  the  trailing  edge.  A final  coat  of  plastic  sealer  was 
applied  to  insure  a minimal  effect  of  line  ink  "thick- 
ness" on  the  uniform  surface  roughness. 

With  an  8 inch  chord  and  12  inch  semispan  (half 
aspect-ratio  of  1.5),,  the  wing  was  attached  to  the 
tunnel  wall  with  a large  endplate  and  protruded  to  the 
center  of  the  test  section.  The  angle-of-attack  was 
varied  to  within  ±1/4  degree  by  a rotation  mechanism 
outside  the  tunnel. 

XX . 3 . Measurement  of  Fluid  Velocities 

A constant  temperature  hot-wire  anemometer  was 
employed  for  the  measurement  of  both  mean  and  turbulent 
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fluid  velocities.  Turbulence  properties  were  measured 

34 

with  an  anemometer  based  on  a design  by  Preymuth  which 

is  detailed  in  Figure  5.  Utilizing  only  three  active 

elements  in  the  sensing  circuit,  this  instrument  exhibits 

low  noise  characteristics  (0.02  microvolts/°C  RMS  mea- 

35 

sured  at  the  input) . An  anemometer  designed  by  Bank 
utilizing  an  operational  amplifier  as  the  error  sensing 
element  was  employed  for  the  majority  of  measurements 
of  the  mean  flow  field  (see  Figure  6) . Though  somewhat 
noisier  than  the  transistor  version,  this  circuit's 
d.c.  stability  over  long  time  periods  (.less  than  2 mV 
per  hour  at  constant  speed  equilibrium)  rendered  it 
extremely  advantageous  for  the  measurement  of  the  poten- 
tial flow  field. 

The  frequency  response  of  these  instruments  was 

obtained  with  a standard  square  wave  response  test. 

The  response  was  never  below  17  kHz  (±3  db)  for  the 

hot-wires  employed  in  this  experiment. 

3 6 

King  investigated  the  heat  transfer  relation 
for  heated  cylinders  in  a uniform  stream.  The  applica- 
tion of  his  results  to  hot-wire  anemometry  gives  rise 
to  a functional  relation  between  the  wire  output  voltage 
and  the  free  stream  flow  velocity — 

E2  - E2  = B Un 


(II-l) 
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where  B is  a constant  depending  on  the  fluid  and  wire 
properties . 

For  the  speed  range  employed  in  this  experiment 

the  exponent  "n"  was  found  to  be  very  close  to  0.5  in 

value  (Figure  7) . As  a result,  an  analog  circuit 

utilizing  Motorola  MC1594  multipliers  in  a squaring 

mode  was  employed  to  linearize  the  anemometer  output. 

The  basic  schematic  is  shown  in  Figure  8.  Squaring 

element  number  2 multiplies  the  anemometer  output  by 
2 2 

itself  giving  E . The  constant  Eq  is  subtracted 

electronically  by  lowering  the  tunnel  speed  to  zero  and 

adjusting  potentiometer,  PI.  The  linearized  output  is 

2 2 

then  obtained  by  squaring  (E  -E  ) with  element  number 
4.  Amplifiers  1 and  3 are  employed  for  gain  control 
making  it  useful  over  a wide  velocity  range.  The 
linearizer  frenuency  response  was  measured  flat  (±3  db) 
to  29  kHz. 

The  hot-wire  probes  utilized  silver  coated 
Pt.  - 10%  Rh.  Wollaston  wire  which  was  soft-soldered 
between  the  ends  of  the  needle  supports.  A weak  nitric 
acid  solution  was  used  in  electro-chemically  etching 
the  silver  jacket  thus  exposing  the  effective  sensing 
element.  The  wire  core  was  0.0001  inches  in  diameter, 
and  wires  were  etched  to  a length/diameter  ratio  of 
approximately  300. 
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Since  X-meters  were  the  primary  probe  types 
employed,  adder  and  subtractor  circuits  were  constructed 

utilizing  Fairchild  Type  741  operational  amplifiers. 

37 

An  X-meter  circuit  constructed  by  McMichael  was 
utilized  for  filtering,  averaging  and  integrating  the 
turbulent  (a.c.)  component  of  the  hot-wire  signals. 

The  circuitry  for  these  elements  is  detailed  in  Figures 
9 and  10. 

Turbulence  intensities  were  also  measured  with 
a Bruel  and  Kjaer  Random  Noise  Voltmeter  (Type  2417) . 
Variable  time  constants  from  0.3  seconds  to  100  seconds 
were  available  through  this  instrument  with  a d.c. 
output  proportional  to  the  RMS  value  of  the  input 
signal.  The  instrument  bandwidth  was  2-20,000  Hz. 

Data  was  recorded  from  digital  voltmeters  and 
a Moseley  type  2D-2A  X-Y  plotter.  A Tektronix  type 
502A  oscilloscope  was  used  for  constant  signal  moni- 
toring . 

The  nature  of  the  investigation  demanded  that  the 
velocity  field  be  resolved  into  components.  A tri- 
axial  3-wire  system  was  considered  but  discarded  due  to 
its  limited  spatial  resolution,  complicated  analog 
circuitry,  and  complex  calibration  procedures.  Since 
an  unseparated,  surface  bounded  flow  was  to  be  studied, 
an  X-meter  type  probe  was  constructed  having  a variable 
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head  angle  (see  Figure  11) . Alignment  of  the  probe 
so  that  the  plane  of  the  "X"  was  locally  parallel  to 
the  wing  surface  effectively  provided  the  total  velocity 
vector  in  most  measurement  regions.  Several  standard 
fixed  angle  X-meters  were  also  employed  for  other  flow 
field  measurements.  The  wire  angle  was  maintained 
nominally  at  45°  inclined  to  the  tunnel  free  stream. 

For  a two-dimensional  flow,  or  a locally  two-dimen- 
sional situation,  the  X-meter  equations  are  given  below. 
Refer  to  Figure  12  for  a graphical  illustration  of  the 
variables : 

Mean  Flow — 


E.  - E_  1 

tan  a = • 

E1  + §2  tan  <j> 

o(E,  + 

U - 1 

2 cos  <}> 


(II-2) 


(II-3) 


o(E9  - E,) 

V = _ 

2 sin  <j> 


Turbulence  Properties — 


T _ °2(el  + e2)2 

U 2 

4 cos  (J> 


(II-4 ) 


( 1 1 — 5 ) 


j?  _ °2(e2  ' ei)2 

. . 2 . 

4 sm  4> 

2 ,~2  T. 

a (e2  - e^) 

4 sin  <{>  cos  <J> 


(II-6) 
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where  the  angles  between  the  free  stream  and  inclined 
wires,  <|)^,  are  equal  and  the  two  anemometer  system 
sensitivities  have  been  adjusted  to  the  same  value,  a. 

The  effective  sensing  region  of  the  X-probes 
was  a rectangular  box  of  dimension  0.02  inches  x 0.02 
inches  x 0.005  inches  (volume  = 2.0  x 10-  inches). 

II. 4.  Probe  Calibration 

Rather  than  rely  on  derived  relationships 

between  the  velocity  components  and  linearized  anemometer 

voltages,  a direct  calibration  of  probe  properties  was 

performed  in  the  tunnel  free  stream.  Employing  the 

proportionality  relationships  (equation  (II-3)  , (II-4)) 

for  the  X-probe,'the  instrument  was  first  aligned  so 

that  its  symmetry  axis  coincided  with  the  free  stream 

direction  and  the  tunnel  speed  varied  to  give  the 

linearized  output  for  "U,"  a typical  value  shown  in 

Figure  13.  The  probe  was  then  rotated  through  known 

angles  in  the  free  stream  at  constant  speed.  A typical 

calibration  curve  is  shown  in  Figure  14.  This  curve 

is  found  to  be  independent  of  the  free  stream  velocity. 

The  validity  of  this  static  calibration  technique  was 

3 8 

verified  by  Singh  who  compared  it  with  a dynamic  cali- 
bration procedure.  Agreement  was  found  to  be  within 
1%.  Figure  14  shows  that  the  crossflow  component,  V, 
is  linear  up  to  a value  of  0.5  U allowing  for  the 
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measurement  of  large  crossflow  components.  Periodic 
probe  calibrations  revealed  that  virtually  constant 
characteristics  were  maintained  over  the  life  of  the 
probe. 

Since  the  wire  angle  with  respect  to  the  free 
stream  is  never  exact  (the  design  angle) , two  calibra- 
tions were  made  for  each  probe.  First,  the  probe 
angle  was  found  for  which  the  inclined  wires  made  equal 
angles  with  the  undisturbed  free  stream.  This  reference 
angle  was  then  used  as  the  basis  of  the  first  cali- 
bration. Secondly,  the  probe  was  aligned  so  that  its 
geometric  axis  coincided  with  that  of  the  free  stream 
velocity  vector.  A second  calibration  utilizing  this 
"zero  angle"  was  obtained.  The  effect  of  this  "forced 
calibration"  is  to  shift  the  curve  to  one  side  or  another 
and  create  an  asymmetric  characteristic.  All  measure- 
ments taken  for  this  experiment,  however,  were  well 
inside  the  linear  bounds  for  all  the  probes  used. 

II. 5.  Measurement  Procedure 

Measurements  of  the  velocity  field  were  taken 
in  both  the  wake  region  downstream  of  the  wing  and 
in  close  proximity  to  the  surface.  For  data  obtained 
in  the  former  region  either  a geometric  or  an  aerody- 
namic (e.g.,  vortex  core  center)  reference  could  be 
utilized.  Measurements  near  the  wing  surface  were  based 
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on  a geometric  reference,  usually  the  wingtip,  estab- 
lished with  the  traverse  positioning  system  in  conjunc- 
tion with  the  optical  cathotometer . 

In  all  cases,  the  probe  was  first  positioned  in 
the  undisturbed  free  stream  with  the  wing  at  the  zero 
lift  angle.  The  linearizers  were  then  adjusted  to  give 
the  desirable  output.  The  probe  was  repositioned  at  a 
point  convenient  to  establishing  the  required  reference. 
For  situations  involving  measurements  close  to  the  wing, 
the  distance  between  the  probe  and  the  surface  was 
established  at  the  position  nearest  contact  by  comparing 
the  cathotometer  readings  for  the  surface  line  reference 
and  sensing  element  locations.  Traverses  were  then  made 
in  horizontal  or  vertical  directions  either  parallel  to 
or  away  from  the  wing.  Graphs  relating  velocity  to 
position  were  obtained  with  the  X-Y  recorder.  In  places 
where  the  turbulence  level  was  significant,  the  data 
were  taken  pointwise  by  time  averaging  the  signal. 

In  addition  to  direct  measurements  of  the 
velocity  field,  information  regarding  the  vorticity 
content  of  various  spatial  regions  was  essential. 

Three  techniques  were  employed  to  obtain  the  rotational 
properties  of  the  flow  field: 

1)  In  the  far  wakn  region,  the  isolated  trail- 
ing vortex  becomes  axisyrjimetric , its  circulation  being 
given  by — 
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T(r)  = 2irr  u0  (II-8) 

A probe  traversal  through  the  vortex  center  with  the 
X-plane  perpendicular  to  the  traverse  direction  allows 
for  the  measurement  of  u0  (r) . Hence,  P can  be  calcu- 
lated using  equation  (II-8) . If  the  vortex  is  com- 
pletely rolled  up,  a constant  value  PQ  will  be  reached. 
Symmetry  was  never  completely  achieved  in  this  experiment 
due  to  the  limitation  on  maximum  downstream  measurement 

location  (x/c  = 7.5  chord  lengths  downstream  of  the 

3 8 * 

trailing  edge) . Singh  indicates  complete  rollup 

(i.e.,  cylindrical  symmetry  of  the  crossflow  velocity 

field)  for  his  higher  aspect  ratio  wing  for  x/c>10. 

2)  In  regions-  of  concentrated  vorticity,  a 

mechanical  probe  similar  to  a "paddle  wheel"  was 

utilized  to  obtain  the  local  value  of  the  vorticity. 

2 6 

The  device  is  described  by  Grow.  Difficulties  in 
using  this  device  in  a quantitative  manner  stem  from 
its  rather  large  effective  volume  (i.e.,  it  must 
average  the  vorticity  over  a relatively  large  region 
in  space)  and  its  inherent  mechanical  friction  which 
limited  its  use  in  weakly  rotational  regions.  It 
provided,  however,  a qualitative  description  of  the 
centers  of  concentrated  vorticity  (e.g.,  vortex  cores) 
and  their  location  around  and  downstream  of  the  wing. 
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3)  A third  technique  employs  the  definition  of 
circulation  as  follows: 

I'  = $ u • dl  = jji  • dA  (IX-9) 

c A 

Contour  integral  paths  were  traced  out  in  physical  space 
by  the  probe  traversing  mechanism.  The  restriction  to 
independent  horizontal  or  vertical  movement  restricted 
these  curves  to  being  rectangles  in  the  Trefftz  plane — 
sometimes  referred  to  as  the  "crossflow  plane"  when  the 
contour  envelops  the  wing.  By  orienting  the  probe 
such  that  the  X-plane  is  parallel  to  the  traversal 
direction,  the  quantities  appearing  in  equations  (II-4) 
and  (IX-9)  were  evaluated.  The  restrictions  and  limi- 
tations of  this  technique  will  be  discussed  below. 

In  all  cases,  the  probes  were  returned  to  the 
zero  reference  location  upon  completion  of  the  measure- 
ment to  examine  the  electronic  drift  of  the  instru- 
mentation. In  situations  where  this  characteristic 
exceeded  one  third  of  1%  based  on  full-scale  values, 
the  data  were  discarded  and  remeasured. 

II. 6.  Restrictions  and  Limitations 
of  the  Measurements 

The  techniques  described  in  the  foregoing 
sections  are  not,  in  general.,  valid  for  all  flow  situ- 
ations. Certain  restrictions  on  the  nature  of  the 


flow  field  are  required  to  insure  compatibility  with 

the  experimental  apparatus. 

The  measurements  taken  near  the  wing  surface 

rely  heavily  on  the  assumptions  of  locally  parallel 

and  unseparated  flow  near  the  surface.  The  first 

assumption  is  invalid  in  the  vicinity  of  the  wing  tip 

due  to  the  surface  discontinuity  that  exists  there  and 

to  the  presence  of  large  rotational  gradients  which 
1 

result  from  the  rollup  of  accumulated  tip-shed  vorti- 
city.  As  a result,  direct  velocity  field  measurements 
in  this  region  are  limited  and  are  presented  only  when 
they  can  be  corrected  for  these  effects  (e.g.,  Quali- 
tative nature  of  the  streamlines). 

Though  not  a problem  in  this  experiment,  the 
effect  of  separated  flow  can  induce  large  errors  into 
the  experimental  results.  Large  reverse  flow  veloci- 
ties with  respect  to  the  probe  body  are  not  readily 
detected  due  to  the  presence  of  the  interfering  probe 
afterbody.  Measurements  upstream  of  the  separation  point 
suffer  errors  due  to  the  non-linear  probe  characteristics 
which  are  significant  for  low  values  of  U due  to 
resulting  large  values  of  V/U  (see  Figure  14) . 

An  additional  difficulty  encountered  in  measure- 
ments neat  a soliu  boundary  with  a hot-wire  anemometer 

is  that  of  probe-wall  heat  transfer.  This  phenomenon 

39 

was  investigated  by  Wills  for  a laminar  boundary 


23 


layer  with  an  extension  to  turbulent  boundary  layers 

40 

and  by  Repik  and  Ponomareva  for  a turbulent  boundary 
layer.  The  closest  point  of  approach  with  the  probes 
used  in  this  experiment  was  0.05  inches  from  the  wing 
surface,  giving  a ratio  of  probe-wall  distance  to 
wire  diameter  of  500.  Wills  showed  that  the  required 
velocity  correction  at  this  distance  is  negligible 
(«  1%) . Data  from  the  other  reference  yield  a similar 
result. 

A major  problem  with  the  use  of  an  X-meter  in 
three  dimensional  flow  fields  is  the  error  due  to  an 
out-of -plane  velocity  component.  This  error  is  cal- 
culated to  be  dependent  on  the  inplane  crossflow 
velocity  component,  the  out-of  plane  component,  the 

wire  inclination  angle  and  the  local  convective  or 
41 

mainflow  speed.  For  small  values  of  the  out-of-plane 
component,  W,  the  linearized  relative  error  to  V is 
computed  as — 

Error  ~ i r 5 — ~ - 7-7  (11-10) 

\o/  sin  ^0  - cos  01 V I 

For  the  standard  45  degree  inclination  employed  in  these 
experiments,  this  expression  reduces  further  to — 


Error  ~ 


(W/U)2 


1 - (V/U) 2 


(11-11) 


It  is  apparent  that  the  boundaries  of  the  contours 
employed  to  determine  circulation  can  greatly  affect 
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measurement  accuracy  if  they  pass  through  a region  where 
large  out-of-plane  component  exists,  e.g.,  regions  of 
concentrated  vorticity.  These  regions  may,  of  course, 
be  enclosed  by  the  contour.  All  experimental  results 
have  been  corrected  for  the  out-of-planeness  error. 


CHAPTER  III 


VELOCITY  FIELD  SURROUNDING  THE 
LIFTING  WING 

III.l.  Optimum  Test  Configuration 

In  addition  to  obtaining  the  parametric  varia- 
tion of  certain  integral  wing  properties,  a primary 
objective  of  this  study  was  to  ascertain  the  detailed 
structure  of  the  flow  field  surrounding  the  lift 
generator  for  some  representative  orientation.  The 
selection  of  an  optimum  operating  condition  was  based 
largely  on  two  considerations: 

1)  High  lift  was  desired  to  accentuate . the 
crossflow  velocities  and  therefore  improve  the  accuracy 
of  lift  related  flow  field  measurements  and, 

2)  A minimal  boundary  layer  displacement 
region  was  sought  in  an  attempt  to  simulate  the  "flat 
plate"  or  idealized  thin  wing'  geometry.  This  condi- 
tion is  usually  optimized  througn  drag  minimization. 

The  angle  of  attack  employed  was  therefore  that  corres- 
ponding to  a maximum  lift  to  drag  ratio. 

Rather  than  relying  on  available  high  Reynolds 
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number  test  data,  the  gross  wing  properties  were 
estimated  through  preliminary  measurements  taken  in  the 


wake.  The  lift  coefficient  was  calculated  by  deter- 
mining the  circulation  about  the  isolated  free  tip 
vortex  downstream  of  the  wing  and  applying  the  approxi- 
mate relation  from  two-dimensional  potential  airfoil 
theory : 

2 r 

CL  = g-f  (III-l) 

o 

Circulation  measurements  were  made  at  7.5  chord  lengths 

downstream  of  the  trailing  edge  using  the  method  out- 
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lined  by  Singh,  Chapter  II.  The  resultant  lift  curve 
is  given  in  Figure  15. 

The  drag  was  approximated  by  evaluating  the 
axial  velocity  defeat  in  the  wake.  Through  the  measure- 
ment of  this  distribution  across  the  wake  (z-direction) , 
the  drag  coefficient  per  unit  span  is  given  approximately 
as — 

2/*Tu  - U ) dz 

r,  _ O X 

D GTc  (I I 1-2) 

The  drag  variation  with  angle  of  attack  is  shown  in 
Figure  16.  This  data  is  based  on  an  average  of  four 
spanwise  locations.  A comparison  of  data-  obtained  at 
several  chordwise  locations  indicated  an  error  of  less 
than  5%. 

The  resultant  expression  for  the  approximate 
lift-to-drag  ratio  is  tihen — 
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The  measured  characteristic  curve  is  given  in  Figure  17. 
This  result  served  as  the  basis  for  the  selection  of  a 
primary  test  angle  of  attack  of  4 degrees.  This  low 
value  allows  one  to  apply  the  small  angle  approximations 
(sin  a *»  a,  cos  a * 1)  in  relating  the  wing-fixed  (or 
"body-fixed")  and  wind  axis  coordinate  systems.  The 
main  flow  axis  here  approximately  coincides  with  the 
airfoil  symmetry  axis. 

The  technique  used  for  obtaining  the  lift 
characteristic  suffers  from  several  deficiencies  in 
this  experimental  situation.  (These  discrepancies  will 
be  discussed  at  length  in  Chapter  V.)  The  method  was, 
however,  sufficiently  accurate  to  reveal  the  shape  of 
the  L/D  characteristic  and  provide  an  adequate  choice 
in  the  selection  of  the  test  orientation. 

A free  stream  tunnel  velocity  of  60  feet  per 
second  was  selected  as  the  maximum  useable  speed  for 
which  probe  vibration  effects  were  negligible.  The 
corresponding  Reynolds  number  based  on  wing  chord  was 
247,000. 

III. 2.  Flow  Regions  in  the 

Crossflow  Plane 

An  examination  of  the  behavior  of  the  flow  field 
in  any  crossflow  (Trefii.:)  plane  intersecting  the  wing 
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surface  reveals  a general  structure  which  is  qualita- 
tively independent  of  the  chordwise  location.  A sketch 
detailing  this  structure  is  given  in  Figure  18.  The 
flow  is  subdivided  into  various  sections  which  exhibit 
similar  properties.  The  relative  physical  sizes  of 
these  regions  will  vary  with  chord  location  and  other 
aerodynamic  parameters.  For  example,  the  thickness  of 
the  inboard  boundary  layer  (Region  2 and  inboard  por- 
tions of  Region  3)  will  vary  with  chordwise  location 
due  to  the  pressure  gradient  (3p/3x)  which  is,  in  turn, 
a function  of  wing  orientation.  The  wavy  lines  which 
depict  the  interface  between  regions  are  not  neces- 
sarily rigidly  defined  but  can  be  a function  of  wing 
construction  uniformity  and  subject  to  some  interpre- 
tation as,  for  example,  is  the  thickness  of  a turbulent 
boundary  layer.  The  advantage  of  this  subdivision 
concept  is  that  it  facilitates  the  description  and 
analysis  within  any  given  region.  While  this  particu- 
lar structure  has  resulted  from  an  investigation  of  a 
specific  rectangular  wing  shape,  the  result  can  be 
generalized  to  a broader  class  of  planforms,  aspect 
ratios  and  angles-of-attack  providing  the  flow  is 
unseparated  in  the  inboard  regions.  A detailed  explan- 
ation of  this  structure  follows  in  the  ensuing  sections. 

The  simplest  region  from  an  analytical  stand- 
point is  that  of  two-dimensional  ("airfoil")  potential 
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flow  inboard  of  the  wing  tip  and  outside  the  boundary 
layers  on  both  surfaces  (Region  1) . It  can  be  located 
about  most  constant  chord,  unswept  wings  of  high  aspect 
ratio  but  more  frequently  is  realized  only  in  the  wing 
symmetry  plane  as  is  thd  case  in  this  experiment.  An 
adjacent  two-dimensional  boundary  layer  (Region  2)  can 
also  be  present.  This  airfoil-like  flow  field  coincides 
with  the  flattening  of  the  spanwise  lift  distribution 
near  the  wing  centerline. 

Observations  nearer  the  tip  indicate  that  a 
spanwise  velocity  component  is  introduced  resulting 
from  the  spanwise  pressure  gradient.  This  component 
becomes  increasingly  more  significant  as  the  tip  is 
approached . 

III. 3.  Characteristics  of  the  Three- 

Dimensional  Boundary  Layer 

The  three-dimensional  wing  boundary  layer 
(Region  3 in  Figure  18)  is  a complex  structure  whose 
velocity  field  is  a function  of  all  spatial  variables. 

It  is  this  region  that  is  solely  responsible  for  the 
generation  of  vortici'ty  in  the  main  flow  direction. 

Measurements  of  the  mean  velocity  components 
and  turbulent  intensities  were  made  at  various  locations 
adjacent  to  the  wing  to  establish  the  extent  of  the 
displacement  region.  The  results  of  several  represen- 
tative measurements  are  graphically  depicted  in 
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Figures  19-27.  In  these  Figures  the  coordinate 
reference  for  the  surface  normal  direction  (z/c  = o) 
is  the  local  surface  location.  The  data  contained  in 
these  illustrations  refer  to  upper  surface  conditions 
only.  The  lower  wing  surface  exhibited  an  extremely 
thin  boundary  layer  (less  than  0.03  inches)  preventing 
any  detailed  measurements  in  its  interior. 

Figure  19  shows  the  variation  of  mainflow 
velocity  ,(0.  ) profiles  with  spanwise  location.  Figure 
23  shows  the  variation  of  the  same  parameter  with  chord- 
wise  position..  The  local  upper  surface  boundary  layer 
thickness  is  seep,  to  be  a strong  function  of  chordwise 
location  but  almost  independent  of  spanwise  position 
except  in  the  immediate  vicinity  of  the  tip.  The  weak 
spanwise  variation  is  probably  due  to  irregularities 
in  wing  thickness  and  surface  roughness  related  to 
construction  tqlerances.  An  examination  of  the  span- 
wise  variation  of  turbulence  intensity  (Figure  27) 
indicates  that  this  may  be  the  case. 

The  layer  thickness  varied  from  immeasurably 
thin  (less  than  .03  inches,  z/c  = .0083)  upstream  of 
the  midwing  point  to  over  0.30  inches  near  the  trailing 
edge.  Measurements  of  Uz/Ux  near  the  surface  verify  that 
the  measured  value  of  the  downwa*  *.  coincides  well  with 
the  calculated  value  of  the  local  surface  slope.  A 
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representative  example  of  these  data  are  given  in 
Figure  28  at  a chordwise  location  near  the  trailing 
edge.  This  provides  some  evidence  as  to  the  integrity 
of  the  mainflow  measurements  near  the  surface  and 
therefore,  the  absence  of  separation.  The  absence 
of  separation  war-  confirmed  by  measurements  over  the, 
angle  of  attack  range  of  0-0  degrees.  Further  discus- 
sion of  this  evidence  will  be  found  in  Chapter  V. 

The  boundary  layer  characteristics  are  observed 
to  be  almost  independent  of  span  location  over  a large 
portion  of  the  inboard  region,  as  noted  in  Figures  22 
and  27.  The  tip  region,  however,  is  characterized  by 
gross  spatial  non-uniformity  in  all  variables.  The 
crossflow  profiles  detailed  in  Figures  20,  21  and  24 
dramatically  illustrate  this  variation.  These  charac- 
teristics will  be  discussed  fully  in  the  next  section. 

The  presence  of  a highly  frequency  selective 
waveform  in  the  fluctuating  velocity  trace  was  observed 
at  the  outer  edges  of  the  boundary  layer  over  the  inboard 
region  near  the  aft  portion  of  the  wing  (x/c>0.70) . 

These  "Tollmien-Schlicting"*  waves  are  associated  with 
the  process  of  transition  to  a fully  turbulent  flow. 


*These  waves  were  first  suggested  by  Rayleigh 
examined  extensively  analytically  by  Tollmien44  and 
Schlicting^5  and  verified  experimentally  by  Schubauer 
and  Skramstad4^  for  the  case  of  flow  over  a flat  plate. 
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Experimental  evidence  of  this  phenomenon  is  displayed 
in  the  correlation  function,  Figure  26.  These  waves 
were  observed  through  the  entire  angle  of  attack  range 
of  unseparated  conditions  (0-8  degrees) , though  the 
amplitude  of  the  waveform  was  greatest  at  the  lower 
angles . 

As  indicated  above,  the  lower  wing  surface 
exhibited  a much  thinner  displacement  region  than  upper 
surface.  Only  in  the  immediate  vicinity  o."  the  cusped 
trailing  edge  did  a suddenly  high  turbulence  level  and 
shallow  mainflow  velocity  gradient  indicate  possible 
local  separation  (or,  incipient  transition) . The  char- 
acteristics of  the  outer  flow  did  not  seem  to  be 
affected  by  this  situation. 

A diagram  depicting  the  growth  of  the  displace- 
ment region  over  the  inboard  portion  of  the  wing  is 
given  in  Figure  25.  This  scale  illustration  shows  the 
chordwise  variation  of  boundary  layer  parameters  for 
a = 4°. 

III. 4.  Inboard  Potential  Flow  Region 

A region  of  considerable  interest  is  that  of 
three-dimensional  potential  flow  just  outside  the 
boundary  layer  inboard  of  the  tip  (Region  4 in  Figure 
18) . It  is  this  sector  that  will  be  involved  in  the 
computation  of  wing  bound  vorticity  in  Chapter  IV. 
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Continuous  spanwise  traverses  were  made  at 
various  distances  from  the  upper  and  lower  surfaces  to 
provide  a complete-  quantitative  description  of  this 
portion  of  the  flow  field.  Both  the  mainflow  (x)  and 
crossflow  (y)  components  were  measured.  Data  were 
obtained  in  5%  chord  intervals  from  the  leading  edge 
to  the  30%  chord  point  and  in  10%  chord  increments  from 
there  aft  to  the  trailing  edge.  Figures  20,  21,  24  and 
29-35  contain  a series  of  various  crossflow  velocity 
profiles  adjacent  to  both  wing  surfaces. 

Examination  of  the  crossflow  variation  reveals 
that  it  is  apparently  not  dependent  on  the  normal  (z) 
coordinate  in  this  region.  This  important  simplifica- 
tion is  more  readily  understood  by  examining  the  equation 
for  the  x-vorticity  component: 
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The  left  side  of  this  expression  is  zero  since  only 
potential  flow  is  considered.  Since,  for  the  untwisted, 
noh-porous  wing  with  a thin  uniform  boundary  layer  the 
streamlines  approximately  follow  the  surface  contour, 
then 
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and  the  above  equation  reduces  to — 
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This  equation  is  verified  by  the  existing  experimental 
evidence.  The  approximation  breaks  down  several  inches 
away  from  the  surface  due  to  the  influence  of  other 
elements  of  the  flow  field.  This  simplification  does, 
however,  allow  a wide  margin  for  error  in  the  measure- 
ment location  employed  for  calculating  the  bound 
vorticity,  as  will  be  noted  in  the  next  chapter. 

The  spanwise  crossflow  velocity  profiles, 

Figures  29-35,  represent  a complete  survey  of  the 
velocity  field  in  the  potential  region  adjacent  to  the 
inboard  upper  and  lower  surfaces  of  the  wing.  Each  data 
point  has  resulted  from  a spatial  average  (in  the  z- 
direction)  in  velocity  for  the  region  for  which  equa- 
tion (III-5)  is  valid.  Error  bars  have  been  added  for 
those  locations  where  the  range  of  values  exceeds  that 
represented  by  the  symbols. 

One  disturbing  feature  that  stands  out  in  these 
data  is  the  reversal  of  crossflow  near  the  trailing 
edge  on  the  upper  surface  and  the  leading  edge  on  the 
lower  surface.  Though  it  occurs  well  inboard  of  the 
tip  and  is  small  in  magnitude,  this  reversed  secondary 
flow  seems  to  oppose  the  expected  spanwise  pressure 
gradient.  Two  explanations  for  the  existence  of  this 
phenomenon  seem  possible: 
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(a)  Local  separation  within  the  boundary  layer 
is  affecting  the  outer  flow  in  the  cross- 
flow component,  and 

(b)  The  tunnel  wall  boundary  layer  is  imposing 
some  effect  on  the  outer  flow  near  the  wing- 
wall  attachment  point. 

Explanation  (a)  does  not  seem  reasonable  in  light  of 
the  results  of  downwash  and  mainflow  measurements  in 
this  region  (see,  for  example.  Figure  28).  Additional 
measurements  (Chapter  V)  which  relate  the  crossflow 
distribution  with  angle-of-attack  also  invalidate  this 
explanation.  The  effect  of  the  tunnel  wall  boundary 
layer,  however,  is  significant  and  may  well  be  respon- 
sible for  this  situation.  Significant  turbulence  (well 
above  the  free  stream  level)  was  detected  as  far  as 
two  inches  distant  from  the  wall  surface.  This  is 
not  surprising  when  it  is  realized  that  the  test  wing 
is  mounted  some  5 feet  downstream  of  the  test  section 
inlet.  The  net  effect  of  this  "thick"  tunnel  boundary 
layer  is  to  reduce  the  effective  aspect  ratio  of  the 
wing.  This  is  done,  however,  in  a rather  non-uniform 
way  as  will  be  seen  in  illustrations  presented  later. 

III. 5.  Flow  Field  Near  the 
Vi'ing  Tip 

The  flow  immediately  adjacent  to  the  wing  tip 
is  characterized  by  the-  presence  of  steep  velocity 
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gradients  over  a small  spatial  region  in  addition  to 

i 

an  irregular  geometric  structure.  Unlike  its  regular  1 

behavior  in  the  inboard  region,  the  crossflow  near  the 
tip  is  observed  to  be  highly  dependent  on  all  spatial 
coordinates  as  can  be  seen  in  Figure  36.  This  area, 
represented  by  Region  5 of  Figure  18,  is  one  character- 
ized by  concentrated  vorticity  which  has  accumulated 
due  to  shedding  upstream  in  the  tip  region.  The 
helical  flow  pattern  so  commonly  associated  with  the 
far-field  isolated  vortices  is  observed  to  develop  well 
upstream  of  the  trailing  edge  as  can  be  seen  in  the 
streak  photograph  in  Figure  37 . Further  examination  of 
qualitative  data  shows  this  "lump"  of  vorticity  to 
adhere  to  both  the  sides  and  upper  surrace  of  the 
tip  region. 

The  actual  identity  of  this  region  is  doubt- 
lessly related  to  the  wing  tip  shape  in  addition  to 
other  wing  parameters.  The  measurements  described  hei~e 
concern  only  a flat  tip  whose  plane  is  perpendicular  to 
the  y-coordinate . Since  the  tip  is  usually  the  source  ■ ' 

of  a singularity  in  theoretical  or  numerical  analyses, 
the  physical  characteristics  of  this  portion  of  the 
flow  field  are  significant. 

The  approximate  shape  of  this  rollup  region  in 
various  Tref ftz-plane  stations  along  the  wing  has  been 
reconstructed  using  detailed  X-wire  data.  Horizontal 
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and  vertical  traverses  at  finely  spaced  intervals  were 
used  to  construct  crossflow  plane  velocity  vectors, 
local  maxima  and  minima, and  crossflow  stagnation  points, 
lines  and  regions.  These  reconstructions  are  shown  in 
Figures  38-43. 

The  swirling  secondary  flow  is  seen  to  orig- 
inate close  to  the  leading  edge  (first  observed  at 
x/c  = 0.05)  and  progress  aft  along  the  side  of  the 
wing.  The  small  scale  and  radius  of  curvature  are 

evident  from  the  sketches.  Its  growth  is  restricted  to 

» 

the  side  portion  of  the  tip  until  the  60%  chord  point 
where  it  apparently  begins  to  "roll  over"  onto  the 
upper  surface.  This  repositioning  movement  is  probably 
caused  by  the  general  influence  of  the  pressure  gradient 
around  the  tip.  By  the  time  the  structure  has  reached 
the  vicinity  of  the  trailing  edge  it  is  positioned 
almost  entirely  above  the  upper  surface.  The  complex 
structure  of  the  secondary  flow  at  various  intermediate 
chordwise  locations  is  caused  by  the  local  wing  tip 
geometry. 

Figure  44  illustrates  more  effectively  the 
regions  where  "reversed"  crossflow  exists,  i.e.,  regions 
where  the  local  crossflow  direction  opposes  that  of  the 
nominal  outer  flow  pressure  gradient.  A line  repre- 
senting locations  of  zero  upwash  (i.e.,  z-velocity 
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component)  is  observed  to  lie  almost  parallel  with  the 
free  stream  extending  rearward  from  the  leading  edge 
stagnation  point.  Measurements  of  its  location  are 
shown  more  accurately  in  Figure  45  which  is  a picto- 
gram  showing  the  locations  of  significant  upwash 
parameters  adjacent  to  the  side  of  the  tip.  As  noted 
in  these  Figures,  this  "stagnation  line"  interesects 
the  upper  surface  location  at  the  approximate  point 
of  rollover. 

It  would  seem  useful  to  be  able  to  ascertain 


whether  vorticity  shed  from  the  lower  wing  surface  could 
roll-over  onto  the  upper  surface.  To  estimate  this 
possibility  consider  a hypothetical  slab,  thickness 
h,  as  shown  in  Figure  46.  Assume  the  pressure  to  be 
uniform  adjacent  to  both  surfaces  having  values  Pu  on 
the  upper  surface  and  P-j.  on  the  lower  surface.  A 
comparison  of  the  characteristic  times  associated-  with 
the  lateral  pressure  gradient  and  the  main  convective 
flow  should  provide  an  appropriate  parameter  for  the 
estimate.  That  is,  consider-- 
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where 

Tlateral — required  for  a fluid  "packet"  to  traverse 
laterally  from  the  lower  to  the  upper  surface  (see 
Figure  46) 
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t . . — time  associated  with  flow  in  the  main- 

convectxve 

stream  direction. 

The  latter  quantity  can  be  estimated  immediately  as — 


Tconvective  Uo  - 

To  calculate  the  former  term  one  requires  that  the  flow 
proceed  around  the  "tip"  in  a regular  fashion.  In  fact, 
assume  that  it  follows  a semi-circular  trajectory  as 
viewed  in  the  Trefftz  plane  (or  in  the  absence  of  any 
main  convective  stream) . With  the  radial  pressure 
gradient  balancing  centripetal  acceleration  to  allow 
this  situation  to  occur,  the  momentum  equation  in  the 
0-direction  can  be  written — 

= P ^ (III-8 ) 

The  average  value  of  the  lateral  flow  velocity  over  the 
assumed  trajectory  can  be  expressed  as — 
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The  radial  pressure  gradient  can  be  estimated — 
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Then — 
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and,  further-- 
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The  presssure  differential  can  be  expressed  in  terms  of 
a lift  coefficient  as  follows: 


Ap 


(III-13) 


Hence, 


The  value  for  x must  be  selected  and  in  this  case 
should  correspond  to  the  length  of  the  slab,  c.  If 
the  resulting  value  of  the  parameter  is  much  greater 
than  unity  rollover  will  not  occur  onto  the  upper  sur- 
face. If  the  value  is  less  than  one  it  will  occur 
upstream  of  the  "trailing  edge."  If  the  value  is  near 
unity  it  should  occur  in  the  vicinity  of  the  trailing 
edge. 

The  case  of  an  actual  wing  is  more  complex, 
however,  and  is  complicated  by  a spatially  variable 
surface  pressure  distribution  and  an  airfoil  shape 
quite  different  from  that  of  the  hypothetical  slab. 

The  geometrical  wing  orientation  also  becomes  an 
important  factor.  Nonetheless,  if  known  experimental 
parameters  are  substituted  into  the  above  equation  for 
the  experimental  configuration  (data  at  a = 4°,  h/x  = 
maximum  thickness  ratio)  it  is  found  that 


N = 0.57, 
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a value  that  might  be  expected  in  view  of  the  experi- 
mental observations. 

To  account  for  airfoil  geometry  one  would  have 
to  determine  the  appropriate  values  for  h,  x and  A p. 

As  might  be  expected  the  flow  surrounding  this 
region  of  concentrated  vorticity  is  dominated  by  its 
presence  (induced  velocity  via  Biot-Savart  Law) . The 
streamlines  in  the  outer  region  exhibit  a much  more 
well-organized  spiral  character  approaching  a cylin- 
drical geometry  as  viewed  in  the  crossflow  plane 
(Region  6 of  Figure  18). 

The  area  denoted  by  Region  7 in  that  Figure  is 
one  dominated  equally  by  the  concentrated  tip  vortex 
structure  and  the  influence  of  the  inboard  surface 
boundary. 

III. 6.  Streamline  Pattern 

Employing  the  crossflow  and  mainflow  data 
described  above, the  streamline  pattern  adjacent  to  both 
the  upper  and  lower  wing  surfaces  was  reconstructed. 
These  diagrams,  assembled  from  planform  velocity 
vectors  are  shown  in  Figures  47  and  48. 

Both  diagrams  appear  to  indicate  shear  in 
certain  locations.  These  are  most  obvious  near  the 
tip  on  both  surfaces  and  near  the  trailing  edge  on  the 
upper  surface.  No  actual  streamlines  are  drawn  in  the 
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immediate  wein-ihy  of  the  tip  in  Figure  47  due  to  the 
swirling  non.f.planar  flow  pattern  which  exists  there. 

e'  the  streamlines  actually  enter 
the' .displacement  .region  the  apparent  "slippage"  actually 
represents  shear ; or  vorticity  production.  In  regions 
where  the  streamlines  lie  wholly  outside  the  boundary 
layer,  however,  the  apparent  shearing  is  accompanied 
by  a variation  of  mainflow  and  downwash  in  the  normal 
coordinate  direction  yielding  an  irrotational  gradient. 

In  contrast  to  the  violent  variations  observed 
in  the  tip.  region,  the- inboard  flow  exhibits  very  small 
values  for  the  crossflow  angle.  Streak  photographs 
taken  of  this  flow  sector  (Figure  49)  indicate  almost 
undetectably  small  angles.  The  sensitive  quantitative 
nature; of  the  Hot-wire  anemometer  does,  however,  detect 
the  distrabuhioft  accurately  even  at  these  small  angles. 
Measurements,  at  one  inch  in  from  the  tip  (y/b  = .083) 
reveal  a maximum'  crossflow  angle  of  only  6.4  degrees 
(at  x/c  = 0>.  80)  with  the  value  dropping  below  one  degree 
at  y/b  * 0> 4 1 for  the  same  chordwise  station. 
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CHAPTER  IV 


VORTICITY  FIELD  SURROUNDING 
THE  LIFTING  WING 

IV . 1 . The  Concept  of  .-.ound 
Vorticity 

The  subject  of  bound  vorticity  was  discussed 
briefly  in  the  introductory  chapter  (pp.  5-6) . The 
viability  of  this  concept  is  not  as  unphysical  as  it 
might  first  appear.  If  one  regards  the  vorticity  in 
the  wing  boundary  layer  as  being  attached  or  "bound" 
to  the  surface,  the  direct  physical  connection  is 
apparent.  The  points  at  which  the  bound  vorticity  and 
free  vorticity  from  the  trailing  sheet  meet  must  then 
be  a region  of  vortex  shedding.  In  lifting  line  theory 
this  occurs  everywhere  immediately  downstream  of  and 
perpendicular  to  the  lifting  line  itself.  In  actual 
flows  the  shedding  will  occur  at  any  point  where  the 
boundary  layer  separates,  e.g.,  the  trailing  edge,  the 
wing  tip,  and  any  location  on  the  wing  where  a separ- 
ation line  might  exist. 

The  replacement  of  a complex  wing  having  a 
thick  airfoil  cross-section  with  a lifting  line  or 
thin  vortex  sheet  is  found  to  yield  integral  results 
(overall  wing,  characteristics)  in  good  agreement  with 
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experiment.  Only  recently  have  complex  numerical 
techniques  provided  suitable  results  with  regard  to  the 
flow  cield  details  with  discrepancies,  however,  still 
existing  depending  on  the  mathematical  formulation 
employed.  In  this  chapter,  the  detailed  vorticity 
field  around  the  lifting  wing  will  be  examined  to  ascer- 
tain the  significance  of  various  portions  of  the  flow 
field. 

IV . 2 . Measurements  of  Spatial 
Vorticity  Content 

Chapter  II  (pp.  20-24)  alluded  to  the  use  of  a 
hot-wire  anemometry  probe  to  measure  circulation  via 
the  definition  of  circulation.  A more  complete  discus-  ' 
sion  of  this  technique  is  found  in  Reference  41.  While 
this  reference  discusses  measurements  in  unbounded  flows, 
the  technique  can  be  employed  in  the  vicinity  of  solid 
boundaries  as  well. 

Consider  a contour,  c,  in  close  proximity  to 
the  cross-section  of  a lifting  wing  as  shown  in  Figure 
50.  Furthermore,  let  the  y-direction  segments  of  the 
contour  lie  completely  outside  the  upper  and  lower 
surface  boundary  layers.  The  circulation  about  region 
A can  be  expressed — 


r = $ 0*d!  = / [U  (y+Ay,z*)  - U (y,z')]  dz' 

A -h/2  Z z • 
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,y+  Ay 

+ / [U  (yS-h/2)  - U (y*,h/2))dy'  (IV-1) 

y y y 

From  Stoke 's  theorem  this  integral  is  also  related  to 
the  vorticity  as  follows: 

A ^ h/2  y+Ay 

pa  = K K‘dA  = f f ZX  dy'  dz’  (IV*‘2) 

-h/2  y 

By  equating  these  expressions  and  taking  the  limit  as 
Ay  approaches  zero,  the  following  expression  is 
obtained — 
h/2 

/ £ dz ' = AU  (y)  (IV-3 ) 

-h/2  X ' y 

where 

AUy(y)  = Uy(y,-h/2)  - Uy(y,h/2) 

Expressed  in  terms  of  an  average  in  the  z-direction — 

lx  (y;n)  = £ AUy{y)  (IV-4) 

Hence,  a direct  t ^lationship  between  the  boundary  layer 
(or  "bound")  vorticity  and  the  surface  crossflow  dif- 
ferential has  been  established.  If  the  wing  and  boundary 
layer  thicknesses  are  negligibly  thin  (h->-o,  a vortex 
sheet),  the  crossflow  differential  defines  the  strength 
of  the  sheet,  h more  complex  derivation  can  be  found 
in  Karamcheti, ^ pp.  530-531.  The  reasons  for  the 
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selection  of  a thin  wing  with  a small  displacement  region 
are  thus  apparent.  A graph  of  AU^  as  a function  will 
then  give  the  shape  of  the  vorticity  distribution  when 
Az  = h is  kept  constant  along  the  span. 

The  utilization  of  this  technique  to  measure 
vorticity  close  to  the  surface  should  then  give  accurate 
results  over  the  inboard  portions  of  the  wing.  In 
this  region  the  streamlines  closely  follow  the  surface 
curvature  and  the  downwash  distribution,  W,  is  found 
to  be  independent  of  the  y-coordinate.  The  outer  flow 
which  exhibits  a crossflow  component  independent  of  z 
(see  Chapter  III)  provides  a suitable  range  of  adjust- 
ment for  the  value  of  h.  The  applicability  of  this 
analysis  is,  so  far  as  the  vorticity  distribution  is 
concerned,  predicated  on  the  assumption  that  the  wing 
plus  boundary  layer  regions  exhibit  a sheet-like 
structure  of  vorticity  with  finite  thickness,  h = 
o(t  + 6U  + fi  ) . 

As  noted  in  Chapter  III  and  in  later  sections, 
the  flow  in  the  vicinity  of  the  tip  does  not  behave  in 
this  manner  and  in  fact  precludes  the  applicability  of 
the  traditional  concept  of  a vortex  sheet.  The  steep 
swirling  gradients  in  velocity  greatly  enlarge  the 
error  in  X-wire  data  due  to  out-of-planeness  of  the 
flow.  A more  expensive  examination  of  the  flow  field 
in  this  region  is  therefore  required. 
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If  one  considers  a contour  extending  from 
y = to  some  value  of  y inboard  of  the  tip  and  to  the 
left  of  the  wing  centerline  (as  shown  in  Figure  51, 
contour  B) , the  circulation  about  this  closed  curve 
can  be  expressed  through  an  extension  of  equation 
(IV-1)  as  follows; 


r(y)  = rtip  + f AUy  (c'y,)  dy’ 
yo 


then 


(IV-5) 


= AUy(c,y)  (IV-6) 

where  c denotes  the  chordwise  location  of  the  contour, 
in  this  case  at  the  trailing  edge.  These  equations  are 
applicable  over  the  inboard  region  of  the  wing  surface. 
Tti  represents  the  overall  tip  circulation  in  that 
region  where  high  rotational  gradients  are  present. 

An  additional  relationship  can  be  established 
by  considering  both  contours  in  Figure  51.  The  one 
labelled  "A"  is  taken  to  completely  encompass  all  of 
the  bound  vorticity  in  a direction  parallel  to  the 
local  section.  The  resultant  circulation  about  this 
contour  will  be  labeled  Similarly,  contour  "B" 

as  explained  above  possesses  circulation  Fn.  Conserva- 

D 

tion  of  vorticity  requires  that  all  of  the  bound  vortex 
lines  passing  through  the  capping  surface  of  "A”  also 
pass  through  "B"  if  it  extends  far  enough  in  the 
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outboard  spanwise  direction  (-y) . Separation  is 
assumed  to  occur  only  at  the  wing  trailing  edge  or 
adjacent  to  the  side  of  the  wing  tip.  Then — 

rA  = rB  (iv-7) 

If  it  is  further  assumed  that  the  wing  is 
uniformly  thin  and  at  a small  angle  of  attack,  the 
following  approximation  can  be  written — 
c 

rA(y)  = I AUx(x',y)dx'  + Istag  (y)  + lyis  (y)  (IV-8) 

where  Igt  represents  the  contribution  to  the  overall 
contour  integral  in  the  vicinity  of  the  leading  edge 
(stagnation  line)  and  Iv^s  relates  the  contribution  due 
to  the  viscous  wake  near  the  trailing  edge.  The  validity 
of  this  expression  as  it  relates  to  lift  will  be  con- 
sidered in  a later  section.  Comparing  equations  (IV-5) , 
(IV-7)  and  (IV-8) , the  following  is  noted — 

Y c 

I‘  . + J AU  (c,y 1 ) dy ' = J AU  (x',y)dx' 

tlP  Vo  Y ° 

+ I . (y)  + I . (y)  (IV-9) 

TJais  establishes  a relationship  between  the  crossflow 

differential  and  mainflow  differential  velocity  fields. 

The  equivalence  of  the  TA  (at  wing  centerline)  and 

(for  the  entire  semi-span)  contours  was  first  experimentally 

4 8 

verified  by  Page  „and  Simmons  in  1925. 


IV . 3 . Velocity/  Vorticity,  Circulation 
and  Pressure  Fields 
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A relationship  between  the  crossflow  velocity 
differential,  "bound"  vorticity,  and  circulation  was 
established  in  the  preceding  section.  These  quantities 
are,  in  turn,  related  to  the  pressure  field  and,  hence, 
the  section  and  overall  lift  coefficients. 

If  the  Kutta-Joukowski  condition  is  imposed  on 
the  velocity  at  the  sharp  trailing  edge,  the  lift  per 
unit  span  is  known  to  be  pUQr,  where  T is  the  circula- 
tion around  a circuit  encompassing  the  wing  chord  and 
boundary  layer  and  parallel  to  the  free  stream.  The 
selection  of  this  contour  in  an  actual  viscous  flow 
is  sensitive  to  a number  of  parameters,  most  importantly 

the  selection  of  the  region  of  intersection  of  the 
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contour  with  the  wake.  Thwaites,  pp.  192-194, 

indicates  that  the  error  i,n  the  selection  of  this 
location  is  negligible  providing  the  wing  lift-drag 
ratio  is  large.  This  connection  between  contour  inte- 
grals in  fully  potential  flow  and  actual  viscous  flow 
was  actually  first  discussed  by  G.  I.  Taylor  in  a note 
to  Reference  50. 

If  this  is  the  case, the  section  lift  coeffi- 
cient can  be  computed  as — 

2 rA(y) 

V 


°1(y)  “ 


(1V-1Q) 
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The  above  is  essentially  the  same  expression 
employed  for  estimating  the  lift  of  the  entire  wing  as 
outlined  in  Chapter  III.  In  that  case  r (y)  was  replaced 
by  The  wing  lift  coefficient,  however,  is  not  given 

by  the  above  expression  but  is  obviously  dependent  on  the 
nature  of  the  T-distribution.  Figure  52  depicts  more 
accurately  the  actual  physical  situation.  As  bound 
vorticity  is  shed  at  the  trailing  edge  the  section  lift 
coefficient  decreases  toward  the  wing  tip. 

The  resulting  value  for  the  parameter  can  then 
be  expressed  as — 

cl  - h = I /b°iay  (IV-11) 

o 

then — 

Id 

cl  - iTH  • E / rA‘y>d*  * rc  W (IV-12) 

o o o 

The  error  involved  in  assuming  a flat-rectangu- 
lar distribution  and  say,  an  elliptic  lift  distribution 
is  (1  - ir/4)  or  about  21%.  The  error  involved  in  the 
selection  of  the  optimum  lift-drag  ratio  angle  of  attack 
using  the  incorrect  lift  estimate  is  discussed  in  the 
Appendix . 

The  expression  for  the  lift  coefficient  can 
then  be  written — 
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cl  = rib  (/  yjax-ay  + I isl  (y')ay 

o o o o 

b 

+ / Ivis  (Y')dy'}  (IV-13) 

o 

The  corresponding  expression  substituting  Tg 
as  given  by  equation  (IV-9)  is  not  properly  posed  since 
the  distribution  of  vorticity  near  the  tip  is  expressed 
as  a single  integral  term, 

A general  expression  involving  the  "B"-type 
contours  can,  however,  be  developed  by  considering  a 
more  general  form  of  these  contours.  Consider  the  two 
general  circuits  shown  in  Figure  53.  These  contours 
do  not  intersect  at  the  trailing  edge  but  at  some  point 
(x,y) . Again,  they  are  assumed  to  encompass  the  entire 
displacement  region.  Again,  the  net  vorticity  passing 
through  capping  surface  "A"  is  equivalent  to  that  of 
"B" . This  is  true  at  least  for  the  case  where  the  lift 
distribution  is  decreasing  monotonically  toward  the  tip 
(see  Figure  52) . Then,  employing  previous  results — 
x 

r(x,y)  = / AUx(x',y)dx'  + lgtag(x,y)  + lAB(x,y) 

= rtip(x,yo}  + f AUy(x,y')dy'  + IAB(x,y)  (IV-14) 
yo 

which  holds  for  points  over  the  inboard  region  on  the 
wing.  The  term,  IAB,  represents  a contribution  to  the 
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contour  inteyral  from  the  downwash  at  point  (x,y) 
where  the  contours  intersect. 

The  lift  coefficient  can  be  established  by 
considering  that  group  of  "B"-contours  which  extend 
all  the  way  to  the  wing  centerline,  y = b.  Realizing 
that  equation  (IV-12)  merely  expresses  the  lift  coef- 
ficient as  an  average  value  of  the  circulation  over 
the  wing  surface,  an  equivalent  expression  can  be 
written — 

cl  = oih?  (x')a*'  (IV-l5) 

o o 


Then — 

Cr  = u T7  if  <x',y«)dy-dx-  + JCr.(x’,y0)dx' 

L uoD  c o yn  y o zip 


+ / IAB(x' ,b)dx' } 


(IV-16) 


Comparing  equations  (IV-12)  and  (IV-15)  , one  finds — 


Cr  = 


'L  Ub  ‘B 
o 


rn(x)  = 


U C A 
o 


(IV-17 ) 


where  the  averages  refer  to  the  spatial  coordinates. 

Since  the  lift  is  simply  the  surface  integral 

of  the  pressure  distribution,  it  is  found  that — 

c b 

/ I’  (x')dx'  = J r,(y’)dy’ 


= pO-  \C  jAp(x',y’)dy'dx' 


o o o 


(IV-18 ) 
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If  two  contours  a distance  Ay  apart  are  con- 
sidered (see  Figure  54),  the  net  lift  on  the  enclosed 
strip  can  be  computed — 
y+Ay 

AL  = J pUQ  r (y')dy'  (IV-19) 

y 

and — 


y+Ay  c 

AL  = f J Ap (x* ,y 1 )dx*dy ' (IV-20) 

y o 

One  can  then  deduce — 

1 c 

VY)  = pir  / Ap(x',y)dx'  = rB(c,y)  (IV-21) 

Similarly,  it  is  also  found  that — 

1 b 

rB(x)  = / Ap{x,y ' )dy 1 (IV-22) 

■ o o 

The  above  relationships  (with  the  exception  of  equation 
(IV-21) ) are  not  strictly  valid  due  to  the  breakdown 
of  the  concept  of  the  sheet  near  the  tip.  They  do, 
however,  provide  a good  approximation  which  relates  the 
velocity  (circulation)  and  pressure  fields  in  an  integral 
fashion. 


Lastly  the  crossflow  at  the  trailing  edge  over 
the  inboard  region  can  be  related  to  .the  pressure  field 
through  equations  (IV-6)  and  (IV-21)  — 


iuy(c,y) 


d 

ay 


c 

J Ap (x 1 ,y)dx' 
o 


(IV-23) 
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IV . 4 . Bound  Vorticity  in  the 
Inboard  Region 

The  surface  crossflow  differential  distribution 
respresenting  the  x-component  of  bound  vorticity  over 
the  inboard  portion  of  the  wing  is  presented  graphically 
in  Figures  55-57 . Spanwise  profiles  are  shown  for  10% 
chord  intervals  along  the  surface: 

Values  of  the  crossflow  velocity  component  have 
been  rendered  non-dimensional  through  division  by  the 
free  stream  velocity.  The  resulting  parameter  can  be 
shown  to  be  independent  of  the  free  stream  value  at  the 
trailing  edge,  i.e., 


AU  (c,y) 


c dciM 


(IV-24 ) 


This  result  follows  from  the  equations  discussed  above. 

The  parameter,  h,  from  equation  (IV-4)  has  been 
adjusted  to  correspond  with  the  maximum  observed  value— 


h = (t  + 6 + 5, ) 

u l max 


(IV-25) 


The  experimentally  determined  value  was  found  to  be  0.72 
inches  or  h/c  = 0.09. 

Examination  of  the  Figures  shows  that  the  net 
value  of  the  crossflow  reverses  polarity  well  inboard 
of  the  tip,  but  not  at  the  tunnel  wall  as  one  would  hope. 
As  indicated  in  the  preceding  chapter  this  is  probably 
the  effect  of  the  tunnel  wall  boundary  layer  on  the 
outer  flow;.  More  e,x'pl  'cdtly,  it  is  the  result  of  the 
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interaction  of  the  tunnel  wall  boundary  layer  with  the 
wall-bound  tip  vortex.  Had  the  wall  boundary  layer  been 
negligibly  thin  one  might  have  expected  this  vortex  to 
be  "suppressed"  with  respect  to  the  external  flow.  In 
this  case,  however,  the  tunnel  wall  boundary  layer 
thickness  is  appreciable  (nearly  two  inches)  so  that 
the  complex  interaction  may  well  influence  the  outer 
f low . 

A more  comprehensive  picture  of  this  portion  of 
the  flow  field  can  be  obtained  from  Figure  58  which  is 
a planform  map  of  the  crossflow  differential  distribu- 
tion. It  can  be  seen  that  the  actual  effect  of  the 
crossflow  direction  reversal  is  to  reduce  the  effective 
aspect  ratio,  though  in  a highly  irregular  manner. 

Experimentally  determined  values  of  the  differ- 
ential can  be  compared  with  theory  through  the  use  of 
the  circulation  function  and  equation  (IV-6)  which  is 
restated  below  in  more  general  form — 

MJy(x,y)  = r (x,y)  ] (IV-26) 

In-  the  case  of  lifting  line  theory,  the  elliptic 
lift  distribution  gives  rise  to  a trailing  edge  value 
of  the  crossflow  differential  which  is — 

ro  (b-y) 

AUy(c'^ elliptic  - (IV-27) 

b\/2by-y2 
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If  the  Schrenk  approximation  is  employed,  the  corres- 
ponding prediction  for  a rectangular  wing  becomes — 


AUy (c,y) 


Schrenk 


1 o (b-y> 

2b\/2by-y2 


(IV-28) 


exactly  half  the  above  value.  Other  values  can  also 
be  predicted  provided  that  the  circulation  function 
is  known. 

These  "predicted"  values  are  compared  with 
measured  data  at  the  trailing  edge  in  Figure  59.  Two 
predictions  are  shown  for  the  Schrenk  approximation. 

The  solid  line  represents  the  value  for  a wing  having 
a twelve  inch  semi-span  which  is  expected  on  the  basis 
of  the  experimental  wing  geometry.  The  other  curve  is 
determined  for  a wing  of  8.5  inch  semi-span  which 
corresponds  to  the  location  where  the  crossflow  differ- 
ential reverses  at  the  trailing  edge.  The  elliptic 
prediction  (dashed  line)  is  determined  only  for  the 
latter  case. 

Actual  values  are  shown  to  lie  between  those 
predicted  by  the  Schrenk  and  elliptic  distributions, 
erring  by  as  much  as  40%.  That  the  experimental  data 
generate  a curve  closely  resembling  that  of  the  predic- 
tors (in  shape)  is  shown  in  Figure  60.  Here  the  cross- 
flow  is  non-dimensionalized  through  division  by  the 
value  at  one-irtch  inboard  of  the  tip  (y/b  = .083). 


! 

( 


1 
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Hence  the  approximate  curve  shape  is  correctly  estab- 
lished. 

The  use  of  the  crossflow  difference  as  a basis 
of  comparison  provides  a more  sensitive  indicator  as  to 
the  detailed  nature  of  the  bound  vorticity  than  the 
more  commonly  employed  circulation  (or  lift)  distribu- 
tion. 


An  additional  question  to  be  addressed  by  these 
results  concerns  the  nature  of  the  circulation  function 
and  the  possibilities  of  a spatially  separable  solution. 
That  is,  can  one  express — 

r (x,y)  = A (x)  B(y)  (IV-29) 

Since  many  analyses  deal  with  the  span  load  distribu- 
tion exc2 isively,  one  might  ask  whether  it  is  reason- 
able to  decouple  the  chordwise  dependence  in  achieving 
such  a solution.  If  this  is,  in  fact,  the  case,  then — 


AU  (x, y)  . , . 
y 1 J _ A (x) 

AUy(c,y)  A (c) 


(IV-30) 


which  is  a function  of  x-only. 

Experimentally  determined  values  of  this  function 
are  shown  in  Figure  61  for  several  different  spanvyise 
locations.  The  parameter  appears  to  be  touch  more 
dependent  on  "x"  than  span  though  its  y-dependence  is 
apparently  not  negligible  on  the  basis  of  these  results. 
Equation  (IV-29)  might  be  more  appropriately  written  as — 
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T (x,y)  = A (x, y)  B (y)  (IV-31) 

IV . 5 . Vorticity  Content  in  the 
Tip  Region 

The  presence  of  large  rotational  gradients  and 
the  complex  distribution  of  concentrated  vorticity  in 
the  tip  region  prohibit  the  application  of  the  simple 
measurement  technique  described  above.  A more  compre- 
hensive description  of  spatial  vorticity  content  was 
therefore  obtained  through  the  use  of  rectangular 
contours  in  the  Trefftz  plane  near  the  wing  tip. 

A series  of  these  contours  forming  a rectangular 
mesh  was  constructed  about  the  tip  and  measurements  were 
made  of  the  crossflow  velocity  component  tangent  to  the 
local  path  direction.  Integration  of  the  velocity  field 
about  the  respective  paths  then  yielded  the  circulation. 

Figures  62-68  show  the  results  of  these  measure- 
ments for  the  test  condition.  Data  were  obtained  in 
chordwise  increments  of  10%  from  the  leading  edge  to 
the  trailing  edge.  The  overall  contour  size  was  deter- 
mined by  the  amount  of  spanwise  distance  influenced  by 
the  "rolled-up"  portion  of  vorticity  as  well  as  the 
spatial  extent  of  that  vorticity.  The  distribution 
and  sizes  of  the  measurement  cells  were  based  on  an 
attempt  to  minimize  the  effects  of  an  out-of -plane  flow 
component  on  the  "X"-wire  probe.  The  streamline 
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reconstruction  technique  and  results  discussed  in 
Chapter  III  were  used  in  making  this  determination. 

A comparison  of  Figures  62-68  and  corresponding 
Figures  38-43  indicates  that  the  regions  of  secondary 
swirl  adjacent  to  th'e  surface  do,  in  fact,  exhibit 
high  concentrations  of  vorticity.  Additional  traverses 
in  this  area  indicated  that  the  outer  edges  of  the  swirl 
zones  roughly  correspond  to  the  boundary  of  the  vorti- 
city bearing  regions  as  one  would  expect.  (A  similar 
correlation  was  established  for  the  inboard  regions  of 
the  wing  with  the  additional  correspondence  between  the 
"vorticity  layer"  and  the  boundary  layer  turbulence 
envelope . ) 

The  production  of  vorticity  (x-component)  is 
observed  to  begin  almost  at  the  leading  edge  and  pre- 
dominantly on  the  lower  wing  surface.  The  concentration 
of  vorticity  there  exceeds  that  adjacent  to  the  opposing 
surface  from  the  leading  edge  downstream  to  the  40% 
chord  point.  From  there  aft  the  opposite  ratio  is 
observed.  Figure  69  depicts  this  ratio  as  a function 
of  ciiordwise  location  for  spanwise  regions  extending 
from  the  tip  to  one  half  and  one  inch  inboard,  respec- 
tively. The  ratio,  T / V . , is  observed  to  be 

independent  of  the  spanwise  location  forward  of  the 
60%  chord  location.  The  deviation  of  the  two  curves 
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downstream  of  this  point  is  probably  due  to  the  "roll- 
over" of  free  vorticity  onto  the  upper  surface  and  its 
interaction  with  the  surface  boundary  layer..  This  inter- 
action is  discussed  more  extensively  in  a later  section. 

The  shedding  of  vorticity  into  the  region  oppo- 
site to  the  side  of  the  tip  is  also  observed  to  occur 
close  to  the  leading  edge.  The  significance  of  the 
vorticity  content  in  this  area  is  shown  in  Figure  70 
which  depicts  the  percentage  of  circulation  associated 
with  this  region  compared  to  the  overall  wing  value  as 
a function  of  the  chordwise  location.  The  amount  is 
observed  to  increase  steadily  from  the  leading  edge  to 
the  40%  chord  point  and  thereafter  remain  constant  to 
the  trailing  edge.  This  observation  does  not  imply  a 
sudden  discontinuance  of  the  shedding  process  but 
rather  that  a kind  of  equilibrium  has  been  established. 
That  is,  the  vorticity  which  is  shed  into  the  side 
region  is  equal  in  amount  to  that  rolling  over  onto 
the  upper  surface  and  out  of  this  region. 

A graphic  illustration  of  the  comparative  amount 
of  shed  vorticity  is  given  in  Figure  71  which  shows  the 
ratio  of.  circulation  in  the  tip  region  to  the  local 
chordwise  value  for  the  entire  w_ing  (T^ (x,o)/rB (x,b) ) 
again  as  a function  of  chordwise  location.  The  resulting 
relation  reflects  a comparison  of  shed  vorticity  to  that 
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being  produced  over  the  surface.  The  curve  is  seen  to 
be  relatively  flat,  decreasing  only  slightly  toward  the 
trailing  edge.  This  indicates  that  the  amount  of 
vorticity  shed  into  the  side  region  nearly  keeps  pace 
with  the  local  production. 

IV . 6 . Effects  of  Tip-Shed  Vorticity 

Vorticity  shed  from  the  lower  wing  surface  into 
the  region  outboard  of  the  tip  has  been  observed  to  roll 
up  as  it  moves  downstream  and  spill  over  onto  the  upper 
surface.  The  movement  of  this  aggregate  of  vorticity 
is  expected  due  to  the  influence  of  the  nominal  pres- 
sure field  near  the  tip.  While  it  lies  adjacent  to  the 
tip  this  vorticity  augments  the  upper  and  lower  surface 
flows  through  a "Biot-Savart"  kind  of  influence.  In 
the  region  where  it  has  rolled  over  onto  the  upper 
surface,  however,  a more  direct  viscous  interaction 
with  the  conventional  boundary  layer  is  responsible 
for  local  deviations  in  the  expected  wing  properties. 

Early  experimentalists  discovered  that  the 

section  lift  coefficient  near  the  tip  of  a wing  could 

deviate  markedly  from  the  theoretically  predicted  value 

depending  on  the  planform  shape  of  the  tip  and  the 

wing's  orientation.  The  case  of  a rectangular  wing 

51 

is  illustrated  in  Robinson,  et  al.,  p.  460,  Figure 
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15  for  various  aspect  ratios.  This  graph  is  reproduced 

here  as  Figure  72  and  shows  a significant  increment  in 

52 

the  lift  near  the  wing  tip.  Pearson  was  the  first 
to  formulate  an  empirical  correction  to  theory  to  account 
for  this  deviation.  When  the  results  of  these  early 
pressure  measurements  are  compared  with  the  flow  "pic- 
ture" established  earlier,  a possible  explanation  for 
the  discrepancy  surfaces.  The  high  velocity  rotating 
flow  above  the  upper  surface  causes  a reduction  in  the 
local  pressure  coefficient  which,  in  turn,  gives  rise 
to  an  incremental  increase  in  the  value  of  the  section 
lift.  This  hypothesis  is  supported  by  measurements  of 
section  pitching,  moment  which  predict  a more  negative 
value  than  expected  in  this  region  and  a corresponding 

chordwise  shift  in  the  center  of  pressure  toward  the 
52 

trailing  edge. 

Actually,  this  result  may  be  viewed  as  some- 
thing of  a limiting  case  of  the  "delta"  wing  which 
experiences  leading  edge  separation.  In  that  case,  the 
resulting  vortex  passes  over  a major  portion  of  the 

upper  surface  and  augments  the  pressure  field  through- 

4 9 

out  (see  Thwaites,  pp.  511-517) . The  comparison 
would  be  much  closer  except  for  the  fact  that  the  tip 
thickness  distribution  controls  the  rollover  location 
in  conjunction  with  the  wing  orientation. 


63 


t 

l 

i 


Unlike  the  case  of  the  delta  wing  with  sepa- 
rated  flow,  the  vorticity  which  rolls  over  at  low  angles 
of  attack  is  not  physically  displaced  from  the  upper 
surface.  Delta  wing  leading  edge  separation  occurs 
predominantly  at  high  angles  of  attack  and  the  result- 
ing vortices  influence  the  surface  conditions  through 
their  induced  field.  In  this  case,  however,  the  close 
proximity  of  the  "free"  vorticity  and  that  of  the 
boundary  layer  results  in  a complex  interaction  which 
is  not  easily  formulated  due  to  the  non-neglectable 
influence  of  viscosity.  It  is  this  interaction  which 
actually  prohibits  the  application  of  the  concept  of 
vortex  sheet  when  modeling  this  region. 

While  many  investigators  have  sought  to  elim- 
inate the  "tip  singularity"  through  mathematical  mani- 
pulation and,  in  some  cases,  judicious  choice  of  the 
. 33 

tip  shape,  few  have  attempted  to  face  the  problem 

on  physical  grounds  and  grapple  with  its  restrictions. 

5 32  53 

Recently,  numerical  analysts  ' ' have  succeeded  in 

verifying  that  rollup  does,  in  fact,  occur  adjacent  to 
the  tip.  Their  results,  however,  do  not  illustrate  the 
"rollover"  phenomenon  since  the  wing  -is  modeled  by  an 
infinitesimally  thin  sheet  of  vorticity  consisting  of  a 
small  and  finite  number  of  bound  filaments.  In  all 
cases,  the  application  of  the  sheet  concept  is  extended 
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over  the  entire  wing  surface  including  the  tip  itself. 

In  these  methods,  regions  of  shedding  are  determined 
a priori  by  the  investigator.  The  results  in  situ- 
ations which  simulate  actual  flows  where  rollover  can 
be  observed  must  be  viewed  with  some  skepticism. 

A primary  question,  then,  concerns  the  predic- 
tion of  shedding  in  the  tip  region  and  the  likelihood 
of  rollover  given  a particular  tip  shape,  specifically 
the  planform  characteristics . If  one  considers  an 
"almost"*  rectangular  wing  having  possible  tip  shapes 
shown  in  Figure  73,  qualitative  statements  can  be  made 
regarding  the  nature  of  the  flow  in  that  region.  Figure 
73a  illustrates -a  possible  situation  where,  under  cer- 
tain conditions,  shedding  will  occur  over  the  entire 
extent  of  the  tip  but  rollover  will  not  occur  due  to 
the  inability  of  the  pressure  gradient  to  force  the 
shed  filaments  over  onto  the  upper  surface,.  In  ^Figure 
73b,  shedding  itself  might  not  exist  at  the  tip  due 
to  the  inability  of  the  spanwise  pressure  gradient  to 
create  a "y"-velocity  component  large  enough.  (At 
high  angles  of  attack,  however,  separation  may  occur 

over  the  tip.)  The  former  case  has  been  numerically 
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verified  by  Rehbach  for  a slightly  swept  configur- 
ation. 

*The  angle,  A,  is  kept  small  enough  so  as  not 
to  effect  the  nominal  rectangular  wing  load  character- 
istics. 


If  the  model  of  the  rectangular  slab  employed 
in  Chapter  III  (Figure  46)  is  recalled,  the  velocity 
around  the  semicircular  secondary  flow  was  found  to 


be — 


U 


0 


Ae 

2p 


1/2 


(IV-32) 


* / 

This  implies  that  the  crossflow  velocities  at  the  tip 
on  the  upper  and  lower  surfaces  are  given  as — 


Uy(h/2)  = - Uy{-h/2)  = J (IV-33) 

Now,  one  considers  the  case  of  a hypothetical,  thin 
slab-wing  at  negligibly  low  angle  of  attack  and  having 
spatially  uniform  upper  and  lower  surface  pressures, 

Pu  and  P-p  respectively.  Then,  since 


L = jAp  dS  = Ap  • S = (Pjl-Pjj)  • S (IV-34 ) 


it  is  readily  observed  that 


Uy  (h/2)  = - Uy(-h/2)  = (IV-35) 

which  is  independent  of  the  slab  thickness  and  overall 
surface  area..  For  shedding  to  occur  (case  "b"),  the 
angle,  A,  must  not  exceed  the  following  value — 


A > tan-1  - tan-1^-^ 


U. 


(IV-36) 


Since  only  small  angles  have  been  considered,  this 
becomes — 


A > 


(IV-37) 
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The  opposite  situation  is  true  for  rollover  (case  "a") . 
The  tip  offset  angle  must  be  greater  than  the  critical 
value  (Equation  IV-37)  for  rollover  to  be  prevented. 

An  additional  restriction  which  surfaces  is  that  of  low 
to  moderate  values  of  angle  of  attack  so  that  the  required 
angle,  A,  does  not  get  too  large  invalidating  the 
original  assumption. 

If  the  restriction  on  a "uniform"  pressure 
field  is  relaxed,  then 


. •,  /Ap(x,o ) 

Critical  " tan  \J— (IV-38> 

If  only  case  "a"  is  considered  (prediction  of  shedding) 
and  rollover  is  neglected, the  pressure  field  at  the 
tip  can  be  written  as  a function  of  the  centerline 
pressure.  Then — 


A 


critical 


tan 


K (b)  Ap (x,b ) 
— 


(IV-39) 


or 


^critical  = tan 


-v 


K(b)  Ap (x,b) 


(IV-40) 


Lastly,  the  case  of  an  actual  airfoil  profile 
is  employed  a 5 shown  in  Figure  74.  In  this  case,  it  is 
readily  shown  that — 


„ / h/o\s.  _ /~HTxj  Ap  (x,o) 
y ' ' \J  2 p t (x) 


(IV-41) 


} 


where  h is  the  local  secondary  flow  (swirl)  diameter 
as  given  in  Chapter  III.  Then, 


67 


A 


critical 


= tan'1  y|&  (lv-42) 


for  the  case  where  rollover  does  not  influence  the 
pressure  field. 

One  can  then  qualitatively  see  that,  for  a more 
complex  planform  like  the  semi-circular  tip  shown  in 
Figure  75,  both  shedding  and  rollover  can  occur  to 
varying  degrees  at  various  locations  along  the  side  of 
the  tip,  depending  on  the  airfoil  shape  and  wing  orien- 
tation. 


IV . 7 . Overall  Wing  Vorticity 
Distribution 

The  quantitative  descriptions  presented  in  the 
preceding  sections  can  be  combined  to  yield  the  total 
distribution  of  vorticity  surrounding  the  wing.  Figure 
76  shows  the  resulting  circulation  profiles  at  various 
chordwise  locations  compiled  from  these  data.  These 
profiles  were  generated  with  equation  (IV-5) . 

All  vorticity  adjacent  to  the  wing  is  included 
in  this  "running  sum"  description  including  the  rolled- 
up  region  discussed  above.  Care  must,  therefore,  be 
exercised  in  attempting  to  relate  these  data  to  the 
surface  pressure  differential  in  the  tip  zone.  The 


absence  of  detail  in  this  region  (in  the  Figure)  is 
intentional  due  to  the  ineffectiveness  of  the  experi- 
mental technique  to  correctly  predict  the  distribu- 
tion there. 

Figure  77  is  a replot  of  the  same  data  but 

referenced  to  the  local  chordwise  maximum,  Y (x,b)  . 

d max 

The  distinctive  similarity  in  these  profiles  indicates 
a corresponding  similarity  in  the  various  chordwise 
values  of  the  spanwise  pressure  gradient.  The  devi- 
ation near  the  leading  edge  reflects  a corresponding 
non-uniformity  in  the  pressure  differential  function 
but  is  also  due  in  part  to  measurement  error.  Large 
curvature  and  radial  streamline  gradients  in  this 
region  made  it  virtually  impossible  to  eliminate  flow 
"out-of-planeness"  errors. 

Figure  70,  referred  to  in  an  earlier  section, 
illustrates  the  comparative  influence  of  various  span- 
wise  regions  on  overall  vorticity  content.  Two  sets 
of  data  points  are  employed  to  display,  overall  wing 
vorticity  development  in  the  chordwise  direction,  both 
reflecting  different  values  of  the  "reference"  circu- 
lation. In  one  case  the  trailing  edge  value,  T , was 
employed  while  in  the  other  the  maximum  observed  value 
(which  occurred  at  x/c  = 0.80)  was  used. 

It  is  of  interest  that  over  half  of  the  vor- 
ticity directed  along  the  x-axis  lies  within  one  inch 
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(y/b  = .083)  of  the  wing.  A significant  portion  also 

appears  as  shed  vorticity  adjacent  to  the  tip.  No 

attempt  has  been  made  to  fit  a curve  through  these 

points  due  to  the  complex  and  unknown  nature  of  the 

b 

chordwise  lift  variation,  P(x)  = J Ap(x,y)dy.  For 

o 

comparison,  however,  the  distribution  predicted  for  a 

two-dimensional  flat  plate  is  also  shown.  A deriva- 
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tion  of  this  distribution  may  be  found  in  Thwaites, 
p.  299. 

The  vorticity  content  of  the  various  spanwise 
segments  exhibits  a different  form  if  referenced  to  a 
"local"  maximum.  Figure  78  illustrates  this  fact  with 
several  representative  graphs  of  rfi  (x,y)/l'B  (c,y)  versus 
chordwise  location  for  several  values  of  span.  The 
vorticity  content  in  various  regions  is  observed  to  be 
nearly  constant  along  the  chord.  This  might  indicate 
that  the  comparable  growth  or  production  of  vorticity 
does  not  vary  significantly  with  chordwise  position 
at  least  for  the  inboard  reaches  of  the  wing. 

The  above  data  can  be  presented  in  a more 
comprehensive  way  and  one  more  familiar  ■ to  proponents 
of  the  various  lifting  wing  theories.  One  considers 
a situation  involving  two  contours,  A^  and  to  the 
left  of  the  wing  centerline  (as  shown  in  Figure  79) 
which  terminate  at  points  and  ' ’ 
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respectively.  Assuming  that  both  contours  yield  the 
same  value  of  the  circulation,  then  the  net  number  of 
equal  strength  vortex  filaments  passing  through  A^ 
must  also  pass  through  . Since  ^A^  and  ^B,,  are 
equivalent,  then 


and  the  net  numbers  of  filaments  intersecting  the 
capping  surfaces  of  and  B2  are  equal.  By  connecting 
surface  points  exhibiting  equal  values  of  rfi(x,y), 
families  of  curves  are  obtained  which  represent  lines 
of  constant  vorticity,  or  the  traditional  "bound  vor- 
tex lines."  It  must  be  emphasized  that  this  description 
is  only  valid  where  the  sheet  concept  can  be  applied, 
i.  e. , interior  to  the  interactive  flow  region  near  the 
tip. 

Figure  80  is  a planform  map  of  the  bound  vortex 
lines  developed  for  the  wing  at  four  degrees  angle  of 
attack.  The  interpretation  of  the  contours  is  given  in  non- 
dimensional  form  for  more  universal  applicability  as  will  be 
discussed  in  the  following  chapter.  An  apparent  dis- 
crepancy indicated  by  the  map  concerns  the  inboard 
movement  of  some  of  the  vortex  lines  near  the  trailing 
edge.  This  tendency,  though  not  governed  by  the  external 
flow  while  the  vorticity  is  "bound,"  seems  to  oppose 
the  motion  at  the  trailing  edge  which  indicates  that  the 
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net  crossfl.ow,  U + U , is  outward  (-y 

^upper  ^lower 

direction) . The  most  probable  explanation  for  this 
difficulty  is  simply  error  in  the  detailed  velocity 
measurements.  Since  the  measurement  of  the  circula- 
tion profiles  is  cumulative,  i.e.,  rB(x,yQ)  is  depen- 
dent on  the  flow  in  the  entire  region  o<y<yQ,  any  small 
errors  in  the  velocity  measurements  are  carried  over 
and  magnified  through  the  computation  of  the  circula- 
tion. If,  in  fact,  these  lines  actually  do  tend 
toward  the  wing  center  line,  possible  explanations 
would  have  to  include  the  lower  surface  separation  at 
the  trailing  edge  and  the  influence  of  the  tunnel  wall 
boundary  layer  as  discussed  above. 

Another  question  which  arises  in  light  of  the 
above  results  is  that  of  an-  equivalent  scheme  employing 
the  "A"-type  contours  alluded  to  earlier.  Figure  81 
depicts  such  a contour  enveloping  a wing  section. 

Measurements  were  made  along  the  surface 

employing  these  contours  (see  equation  (IV-8))  in  an 

attempt  to  reconcile  the  two  methods  of  measurement. 

The  paths  were  selected  so  that  all  bound  vorticity  was 

completely  enveloped  and  the  vertical  'path  segment  near 
« 

the  trailing  edge  (x/c  = 1.05)  was  perpendicular  to  the 
free  stream.  Selection  of  this  wake  portion  of  the 
contour  was  based  on  considerations  discussed  in 
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Thwaites,  pp.  192-194,  and  reviewed  briefly  m 
Section  IV. 2. 

Measurements  of  the  locally  parallel  component 
of  the  mainf low  were  made  in  the  inboard  region  with  a 
representation  sample  of  the  results  shown  in  Figure 
82.  Graphs  of  the  resulting  differential  function 
have  the  form  shown  in  Figure  83.  Visual  inspection 
of  these  curves  reveals  that  a large  fraction  of  the 
vorticity  resides  within  one  inch  of  the  wing  tip. 

The  largest  source  of  error  accompanying  these 
measurements  was  that  in  the  velocity  data  near  the 
leading  edge.  The  presence  of  the  stagnation  point 
and  uncertainties  in  flow  direction  created  difficul- 
ties in  obtaining  reliable  results  in  this  region. 

The  chordwise  variation  in  vorticity  content 
is  obtained  through  the  application  of  equation  (IV-14) . 
The  contribution,  IAB(x,y),  was  not  measured  except  aft 
of  the  trailing  edge  so  all  profiles  must  be  viewed  only 
in  a comparative  light.  Figure  84  illustrates  this  data 
when  nondimensionalized  by  local  spanwise  maximum  values 
Again,  the  predicted  two-dimensional  flat  plate  distribu 
tion  is  presented  for  comparison.  Except  for  the 
curve  representing  y>  - 0.5  inches  (y/b  = .042),  the 
similarity  in  curve;  shape  is  remarkable.  This  evidence 
strongly  supports  the  contentions  regarding  the 


similarity  in  the  pressure  differential  function  and 
vorticity  production  made  above.  Specifically,  the 
observations  of  Figures  70  and  84  indicate  that — 
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rB(x,y) 

rjrxTbT  = p (y)  only 


(IV-44) 


and 


rA(x,y) 

T^TcTyT  = G (x)  only 

Expanding  the  upper  relation  gives 

Ftip(x)  + /y%(x,y')dy«  = F(y)  {Ttip(x) 

^o 

b 

+ / AU  (x,y')dy‘  } 

Vo  • 

Differentiation  with  respect  to  y indicates — 

b 

AU  (x,y)  =F'(y)  J AU  (x,y ' )dy 1 
x y y 

Jo 

The  remaining  integral  is  a function  of  x only, 


(IV-45) 


(IV-46) 


(IV-47 ) 


so 


that — 


AUy(x,y)  = F'(y)‘  R(x)  (IV-48) 

Similar  treatment  of  the  second  relationship  reveals — 
AUx(x,y)  = g'  (x)  H(y)  (IV-49) 

In  both  cases,  the  flow  velocity  differentials  are 
separable  functions  of  the  spatial  coordinates.  It 
must  be  emphasized  that  these  relationships  are  valid 
only  where  curve  similarity  is  observed/  e.g. , 
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inboard  of  the  tip  and  downstream  of  the  leading 
edge. 

These  observations  provide  further  support  for 

the  relationship  posed  by  equation  (IV-29)  over  a large 

portion  of  the  surface.  j 

A comparison  of  the  measurements  applied  to  the 

two  types  of  integrals  is  shown  in  Figure  85.  The 

functions  TA(x,b)  and  rQ(c,y)  are  shown  to  be  in  close 

agreement.  Also  depicted  are  the  distributions  obtained 
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by  DeYoung  and  Harper  — modified  lifting  line  theory, 
the  Schrenk  approximation55  and  an  elliptical  lift 
distribution.  Several  data  points  obtained  in  early 
measurements  of  the  surface  pressure  differential  are 
also  shown. 

The  pressure  data  and  circulation  measurements 
are  seen  to  agree  in  all  but  the  immediate  vicinity 
of  the  tip,  as  expected.  Only  the  Schrenk  approxima- 
tion appears  to  adequately  approximate  the  correct 
behavior  of  the  vorticity  distribution.  Again,  however, 
large  errors  persist  in  the  tip  zone. 


CHAPTER  V 


FLOWFIELD  DOWNSTREAM  OF  THE  WING 

V. 1 . Measurements  in  the  Near  Wake 

In  addition  to  the  investigation  of  the  flow 
field  surrounding  the  wing,  measurements  were  made  in 
the  "near  wake"  region  to  determine  its  physical  pro- 
perties . 

Velocity  field  data  were  obtained  to  a maximum 
of  5 chord  lengths  downstream  of  the  trailing  edge, 
specifically,  at  downstream  locations  of  x/c  = 1.0 
(trailing  edge),  2.0,  4.0  and  6.0.  The  wing  angle  of 
attack  was  also  varied  to  values  of  4,  6 and  8 degrees, 
all  representing  conditions  of  largely  unseparated  flow 
over  the  wing  surface.  The  total  wind  tunnel  test 
section  length  was  responsible  for  limiting  the  maximum 
downstream  measurement  location. 

The  contour  integration  techniques  discussed 
in  preceding  chapters  were  employed  to  obtain  the 
spatial  vorticity  distribution  in  the  wake  region. 
Trefftz  plane  contours  like  that  of  Figure  86  were 
employed  for  the  surveys.  The  z-coordinate  levels 
(horizontal  path  elements)  were  selected  to  entirely 
encompass  the  identifiable  turbulent  portions  of  the 
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Wake,  as  well  as  the  central  vortex  core.  Hence,  all 
wing  generated  vorticity  was  contained  within  the 
contour  bounds. 

Measurements  of  the  wake  shape  (turbulent  edges 
and  mainflow  profiles)  were  also  obtained  in  an  effort 
to  identify  the  physical  shape  of  the  vorticity  package. 

NACA  and  other  experimental  data  were  employed 
for  purposes  of  comparison,  where  possible. 

V. 2 . Trailing  Edge  Vorticity 
Distribution 

Crossflow  velocity  differential  measurements  were 
made  at  the  trailing  edge  (x/c  =1.0)  at  the  three  test 
angles-of -attack.  These  data  were  obtained  outside  of 
the  surface  displacement  regions  but  in  the  quasi-two- 
dimensional  inboard  region  discussed  in  preceding  chap- 
ters. Additional  data  were  obtained  which  verify  the 
distinctive  properties  of  this  region  for  all  values  of 
a employed . 

The  inboard  vorticity  distribution  is  shown  in 
Figure  87.  Again,  non-dimensionalization  of  the  cross- 
flow  parameter  is  accomplished  with  the  undisturbed 
mainflow  velocity.  The  curve  shapes  show  a marked 
similarity  indicating,  similar  local  flow  characteristics. 

This  similarity  is  even  more  graphically  illus- 
trated via  the'  circulation  function,  as  shown  in  Figure 
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88.  Shown  as  distributions  with  respect  to  spanwise 
locations,  these  curves  offer  the  additional  feature 
of  including  the  contribution  of  the  tip  increment  to 
the  overall  vorticity  sum.  The  data  strongly  support 
the  contention  that  the  wing  flow  field  characteristics 
are  independent  of  the  angle-of-attack,  at  least  over 
the  inboard  region.  The  unseparated  nature  of  the 
surface  flows  is  probably  the  key  to  this  relationship. 

Care  must  be  exercised  in  inferring  the  simi- 
larity of  the  resulting  pressure  distributions, however , 
especially  in  the  vicinity  of  the  wingtip.  Though  an 
extensive  survey  of  the  surface  flow  was  not  made  at 
all  values  of  a,  -reasoning  presented  in  earlier  chap- 
ters indicates  that  vortex  "spillover"  onto  the  upper 
surface  should  occur  farther  upstream  with  increasing 
angles  of  attack.  This  would  imply  a higher  value  for 
the  tip  lift  increment  (due  to  interaction  with  the 
surface  boundary  layer)  and  a probable  variation  in 
the  section  center  of  pressure. 

The  inboard  flow  field,  being  independent  of 
a, should  then  be  expected  to  exhibit  the  appropriately 
non-dimensionalized  characteristics  of  the  test  case 
discussed  in  preceding  chapters.  Therefore,  the  flow 
"picture"  presented  in  Figure  80  is  valid  for  other 
values  of  a for  which  the  inboard  flow  is  unseparated. 


78 

Only  the  interactive  tip  zone  characteristics  should 
be  expected  to  exhibit  local  variations. 

V. 3 . Wing  Lift  Characteristics 

Values  of'  the'  overall  wing  lift  coefficient 
were  computed  from  the  data  described  above  and  compared 
with  each  other  as  well  as  other  theoretical  and  experi- 
mental data.  The  calculation  employe.d  equation  (IV-16) 
and  circulation  data  obtained  aft  of  the  trailing  edge 
which  will  be  discussed  below. 

Since  a direct  comparison  with  theory  and 
other  data  was  desired,  the  raw  experimental  data  had 
to  be  corrected  for  various  wind  tunnel  effects.  Theo- 
retical and  empirical  corrections  for  solid  blockage 
(of  both  wing  and  traversing  mechanism)  and  lift  effect 
were  utilized.  A complete  discussion  of  these  and  other 

wind  tunnel  interference  effects  can  be  found  in  Pank- 
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hurst  and  Holder,  pp.  328-364.  The  effect  of  solid 
blockage  was  found  to  be  negligible  at  the  measurement 
conditions  due  largely  to  the  small  projected  cross- 
section  at  these  low  angles  of  attack.  The  correction 
for  lift  effect,  however,  was  more  significant  resulting  in 
an  effective  increase  in  angle  of  attack  of  approximately 
0.5  degrees. 

The  various  curves  of  wing  lift  coefficient 
versus  angle-of-attack  are  displayed  in  Figure  89. 


Linear  curve  fits  were  employed  for  the  experimental 
data  due  to  the  unseparated  nature  of  the  flow  in  this 
case. 

Data  are  presented  which  illustrate  the  error 
incurred  from  the  use  of  the  "1/r"  assumption  for  the 
vortex  behavior  in  the  near  field  (see  section  III.l). 
These  measurements,  corrected  for  the  approximate  load 
distribution,  predict  a lower  value  for  than  any  of 
the  theoretical  schemes  or  the  contour  integral  data. 
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Lifting  surface  theory,  here  due  to  Multhopp, 
is  observed  to  predict  the  lift  curve  slope  with  more, 
accuracy  than  the  lifting  line  theory  as  would  be 
expected . 
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. Numerical  results  obtained  by  Rom,  et  al., 
using  the  iterative  numerical  technique  discussed  in 
Chapter  I come  even  closer  to  the  experimental  data. 
These  results  combine  more  "details"  of  the  flow  field 
behavior  than  previous  investigators  have  employed. 
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An  empirical ■ relationship  attributed  to  Wood, 
seems  to  provide  the  closest  agreement  in  this  case. 
Wood.' s relation  is  baeed  largely  on  an  average  of  a 
large  array  of  empirical  and  experimental  results. 

V.4.  Wing  Wake  Vorticity 
■Distribution 

The  wing  wake.,  whose  detailed  shape  is  discussed 
in  a later  section,  can  be  envisioned  as  consisting  of 
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a sheet-like  structure  extending  beyond  the  trailing 
edge  across  the  entire  span  and  a core  of  concentrated 
vorticity  in  the  outboard  region.  The  formation  of 
this  core  has  been  discussed  in  the  preceding  chapters. 

Of  interest  regarding  the  rollup  process  is  the 
rate  of  transfer  of  vorticity  from  the  inboard  wake  to 
the  core  region,  as  well  as  the  core's  evolution.  The 
former  relationship  was  investigated'  for  the  near  wake. 
The  spanwise  distribution  of  "x"-vorticity  is  displayed 
via  the  circulation  function,  with  integration  from 
y = _oo  toward  the  wing  centerline. 

Variation  of  the  circulation  profiles 
downstream  for  fixed  values  of  a is  illustrated  in 
Figures  90-92.  The  reference  values  of  circulation, 

I’  , obtained  for  the  various  bounding  contours  is  shown 
in  Figure  93.  Agreement  is  within  experimental  expec- 
tations. Vorticity  generated  outboard  of  the  tip  (at 
the  trailing  edge)  is  observed  to  move  toward  the 
"Betz"  value  of  the  spanwise  location.  A significant 
outboard  movement  of  the  "sheet"  vorticity  was  not 
observed,  however.  Figure  94  shows  more  clearly  the 
chordwise  variation  of  circulation  for  several  span- 
wise  locations.  An  increase  in  this  function  with  the 
"x"-cordinate  would  indicate  vorticity  transfer  toward 
the  vorticity  centroid.  No  definitive  statement  about 
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this  trend  can  be  made  over  the  region  of  investigation. 

Only  at  the  highest  angle  of  attack  (8  degrees)  does 

the  expected  behavior  occur. 

Some  comparative  data  obtained  by  Fage  and 
48  . . 

Simmons  using  a similar  measurement  scheme  for  an 
aspect  ratio  6 rectangular  wing  is  displayed  in  Figure 
90.  Their  observations  appear  to  fit  within  the  exper- 
imental scatter  of  the  present  measurements.  This 
result  indicates  that  the  rollup  process  is  only 
weakly  dependent  on  aspect  ratio  for  the  entire  sheet. 

This  conclusion  is  consistent  with  the  results  of 
2 6 

Grow  who  achieved  the  same  qualitative  variation  for 
higher  angles  of  attack. 

Figures  95  and  96  show  the  variation  of  x- 
. vorticity  content  with  angle-of-attack  for  fixed  down- 
stream locations.  T ,e  correlation  between  outboard 
vorticity  transfer  and  increasing  angle-of-attack  appears 
to  be  stronger  than  with  downstream  location. 

V. 5 . Shape  of  the  Near  Wake 

Since  the  "bound"  vorticity  adjacent  to  the 
wing  surface  is  located  within  the  boundary  layer 
(unseparated  displacement  region),  it  is  not  unreason- 
able to  expect  that  the  shed  vorticity  downstream 
should  be  confined  to  the  resulting  wake.  Except  for 
the  properties  of  the  vortex  core  this  region  is 
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identifiable  by  the  usual  characteristics  of  a high 
turbulence  level  and  a perceptible  velocity  defect  in 
the  mainflow  direction.  Both  of  these  properties  were 
studied  to  identify  the  wake  location  corresponding  to 
the  vorticity  measurements. 

Figures  97-99  illustrate  the  location  of  the 
wake  edges  for  the  various  values  of  a.  The  data 
represent  spatial  locations  where  the  turbulence  indi- 

p=2~ 

cator,v  w'  /U  , dropped  below  the  value  0.002,  which 
is  three  times  the  free  stream  value.  Additional  data 
indicating  the  locations  of  the  maximum  x-velocity  (U) 
deficit  are  superimposed  on  these  Figures. 

It  is  observed  that  the  wake  at  the  wing  center 
line  descends  more  rapidly  than  at  any  other  location. 

A marked  curling  of  the  inboard  wake  is  also  apparent. 
The  vortex  core,  on  the  other  hand,  appears  to  extend 
downstream  in  nearly  the  free  stream  direction,  posi- 
tioned slightly  above  the  projected  z-location  of  the 
trailing  edge. 

These  observations  are  in  qualitative  agreement 
w.ith  the  results  of  both  Spreiter  and  Sachs^  and 
Chigier  and  Corsiglia^0  who  experimentally  investigated 
the  wake  shape  in  this  region. 

The  waxe  is  observed  to  be  thickest  in  the 
inboard  region  with  a marked  "necking  down"'  toward  the 
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central  core.  Tight  curling  near  the  core  is  observed 
at  x/c  = 2.0  and  to  a lesser  extent  at  other  down- 
stream locations.  Diffusion  at  the  wake  edges  in 
this  region  tends  to  blur  the  expected  spiral  behavior. 

Care  must  be  exercised  when  employing  these 

results  in  a comparison  with  other  experiments  (actual 

flow  or  numerical)  due  to  the  effects  of  wind  tunnel 

2 

wall  interference.  Hackett  and  Evans  investigated  the 
influence  of  the  proximity  of  the  wall  boundaries  and 
their  geometry  and  noted  that  "Wall  proximity  causes 
tip  vortices  to  sweep  downward  less  and  floor  proximity 
has  an  obvious  effect  on  the  penetration  of  the  center 
of  the  vortex  sheet."  Their  experimental  evidence  sub- 
stantiates this  claim. 

If  one  assumes  that  all  non-core  located 
vorticity  lies  along  a line  representing  the  wake  center 
then  the  preceding  description  of  the  vorticity  vari- 
ation (section  V.4)  along  with  the  geometry  of  this 
line  sufficiently  described  the  wake  rollup.  As  indi- 
cated above,  however,  the  wake  exhibits  a finite 
thickness  which,  in  this  case,  cannot  be  presumed  to 
be  insignificant.  The  z-variation  of  - the  vorticity 
may,  in  fact,  cause  some  modification  in  the  wake 
evolution  and  rate  of  rollup.  This  variation  was  not 
examined  in  the  present  study. 
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V. 6 . Rate  of  Rollup;  Discussion 

The  results  of  the  above  sections  can  be 
combined  to  give  a more  graphic  picture  of  the  rollup. 

The  planform  location  of  vortex  filaments  can 
be  constructed  in  a manner  identical  to  that  described 
in  Chapter  IV  for  the  bound  filaments.  Rad  this  been 
accomplished  one  would  note  only  a slight  spanwise 
displacement  of  the  filaments  over  the  region  of 
interest  at  least  ror  small  angles.  This  indicates, 
perhaps,  a qualitatively  good  reason  for  Prandtl ' s 
placement  of  the  free  vorticity  elements  in  the  lifting 
line  theory. 

One  glaring  conclusion  is  inescapable  from  the 
above  result.  The  wake  rollup  process  is  not  nearly 
as  marked  as  one  would  expect  based  on  the  rapid  evolu- 
tion of  the  core  adjacent  to  the  wing  surface.  Though 
the  near  wake  shape  distorts  and  some  transfer  of 
vorticity  into  the  core  takes  place,  it  is  near  the 
lifting  surface  that  the  most  marked  evolutive  changes 
take  place  probably  due  to  its  ability  to  sustain 
shear. 

A number  of  investigators  have  attempted  to 
predict  the  roll-up  distance  (or  time)  employing  various 
bases. 

Kaden's  early  work  led  to  the  belief  that 
rollup  had  occurred  fu]Jy  at  almost  the  wing  itself. 


i. 


— S5K 
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His  basis  for  the  rollup  time  was  the  coincidence  of 

the  vortex  core  with  the  "Betz"  centroid.  His  analysis, 

however,  was  applied  to  an  open-sided  infinite  vortex 

4 

sheet.  A slightly  modified  analysis  by  Moore  using  a 
spanwise  finite  sheet  predicted  99%  rollup  within 
3.5  chord  lengths  of  the  trailing  edge  using  the  same 
criterion . 

Spreiter  and  Sachs^  slightly  improved  analysis 
employing  the  bound  vortex  and  a semi-inif inite  sheet 
yielded  a value  of  4.45  chord  lengths  for  the  present 
test  case.  Again  the  criterion  was  the  same. 

In  all  of  the  above  cases,  agreement  between  the 
estimates  and  the  current  experiment  is  poor.  Perhaps 
the  data  is  in  error,  but  it  is  more  probable  that  what 
is  at  fault  is  the  rollup  criterion.  It  has  been  shown 
that  a vortex  "center"  does,  in  fact,  coexist  with  a 
sheetlike  structure  of  vorticity  extending  to  the  wing 
centerline. 

A more  realistic  criterion  would  seem  to  be  the 
requirement  of  "1/r"  behavior  of  the  velocity  field, 
though  this  location  is  difficult  to  predict. 

Another  interesting  parameter  is  the  non- 
dimensional  distance  between  the  trailing  edge  and  some 
location,  1*,  where  an  infinitesimally  weak  (no  vor- 
ticity) fluid  element  initially  at  the  wing  centerline 
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has  revolved  around  the  central  core  one  revo'ution.  A 
simplified  analysis  employing  a fully  rolled  up  potential 
vortex  extended  rearward  from  the  tip  at  the  trailing 
edge  gives  rise  to  an  upper  bound  in  this  case.  It  can 
be  shown  that — 

2 2 
4tt  b*  U 

1*  = ri-°  (V-U 

O 

where  b*  is  the  effective  semi-span.  For  an  ellipti- 
cally  loaded  wing  this  can  be  written — 

1*  = llLP*  (V-2) 

CL° 

In  the  present  case,  for  a = 4 degrees,  the  approximate 
value  of  this  parameter  is  1*  = 612,  i.e.,  612  chord 
lengths  downstream  of  the  trailing  edge. 
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CHAPTER  VI 

SUMMARY  AND  CONCLUSIONS 

The  following  list  summarizes  the  major  results 
of  this  investigation  and  its  primary  conclusions: 

1.  The  experimental  technique  employed  for  the 
measurement  of  the  circulation  and  vorticity  distribu- 
tions appeared  to  provide  satisfactory  data  for  the  flow 
conditions  utilized. 

2.  The  flow  field  surrounding  the  wing  could 

be  separated  into  various  spatial  regions  each  exhibiting 
unique  properties.  This  separation  facilitated  the 
local  simplification  of  the  flow  field  from  the  more 
complex  general  three-dimensional  situation. 

3.  The  surface  boundary  layers  were  found  to 
contain  all  of  the  wing  generated  vorticity  except  in 
the  immediate  proximity  of  the  tip.  This  verification 
lends  credence  to  the  application  of  the  vortex  sheet 
as  a replacement  model  for  the  actual  wing  if  one 
envisions  the  surface  boundary  layer  and  wing  thick- 
nesses approaching  zero. 

4.  The  flow  field  near  the  wingtip  is  observed 
to  exhibit  a high  degree  of  rollup  immediately  adjacent 
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to  the  wing  surface  similar  to  that  expected  in  the 
wake.  The  resulting  crossflow  plane  secondary  flow 
is  observed  immediately  downstream  of  the  leading  edge 
in  the  region  opposite  the  side  of  the  tip.  Vortex 
shedding  from  the  lower  wing  surface  is  responsible 
for  this  occurrence. 

5.  This  rolled-up  "free"  vorticity  is  seen 
to  roll  over  onto  the  upper  surface  at  a predictable 
chordwise  location  for  the  simple  geometric  tip  shape 
employed. 

6.  The  behavior  of  the  shed  vorticity  at  the 
tip  can  modify  the  expected  surface  pressure  distribu- 
tion and  resultantly  invalidate  the  applicability  of  the 
vortex  sheet  concept  due  to  direct  viscous  interaction 
with  the  conventional  wing  boundary  layer. 

7.  Investigation  of  the  inboard  bound  vorticity 
distribution  revealed  that  the  bound  circulation  func- 
tion could  be  expressed  approximately  as  the  product 

of  two  separable  functions  of  the  spatial  coordinates, 
x and  y. 

8.  A rolled  up  vortex  with  an  identifiable 
core  was  observed  to  emanate  from  the  upper  surface  at 
the  trailing  edge.  This  region  was  found  to  contain 
over  60%  of  the  total  wing  generated  vorticity  for  the 
configuration  tested. 
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9.  The  trailing  edge  vorticity  distribution 
was  found  to  be  independent  of  the  wing  angle  of  attack 
for  the  unseparated  conditions  studied.  There  is  good 
justification  to  assume  that  the  entire  surface  bound 
vorticity  field  exhibits  similar  characteristics. 

10.  The  wing  lift  coefficient  is  found  to  vary 
nearly  linearly  with  angle-of-attack  as  expected. 
Corrected  experimental  results  are  found  to  be  in  good 
agreement  with  direct  force  measurements  and  various 
theoretical  and  numerical  predictions. 

11.  A comparison  of  near  surface  data  and 
measurements  in  the  near  wake  region  indicate  that 
rollup  occurs  much  more  violently  near  the  solid  surface 
than  in  the  boundary-free  zone  downstream  of  the  wing. 

12.  The  normalized  description  of  the  rollup 
process  in  the  near  wake  region  is  observed  to  be  nearly 
independent  of  both  aspect  ratio  and  angle-of-attack. 

13.  The  isolated  vortex  core  is  observed  to 
move  downstream  in  nearly  the  free  stream  direction 
while  the  inboard  wake  is  displaced  downward  (opposite 
the  lift  direction) . The  largest  displacement  is 
observed  to  occur  near  the  wing  centerline.  A distinc- 
tive curling  motion  is  apparent  as  the  v/ake  evolves. 

14.  The  thickest  portion  of  the  wake  (z- 
direction)  occurs  downstream  of  the  wing  centerline 
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with  a marked  necking  down  observed  near  the  rolled  up 
core.  A complete  examination  of  the  roll  up  charac- 
teristics should  include  variation  in  the  vorticity 
distribution  in  this  direction. 

15.  Various  estimates  of  the  roll  up  rate 
reviewed  in  this  study  are  found  to  err  significantly, 
predicting  a fully  rolled  up  geometry  far  upstream  of 
its  actual  occurrence.  The  "x"-oriented  vorticity  in 
the  inboard  wake  region  is  observed  to  be  ingested  into 
the  core  at  a slow  rate  in  the  near  wake  region. 
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EFFECT  OF  LIFT  CURVE  SLOPE  ON  THE  MAXIMUM 
LIFT-DRAG  ORIENTATION 


The  condition  for  maximum  lift-drag  ratio  is 

simply — 


(L/D) 


0 


For  wings  of  moderate  to  high  aspect  ratio,  the  conven- 
tional relations  are 
CL  = MO 

CD  “ K1  + K2  CL 
= Kj_  + K2  M2  a2 


Then 

K,  1/2 

a = { L_  1 

max  L/D  * ..2  f 

K~M 


A small  error,  AM,  in  the  measurement  of  the  lift  curve 


slope  results  in  an  error  in  a T of — 

max  L/D 


Aa  _ 
max  L/D 


AM 
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FIGURES 


Fig. 


Behavior  of  Fluid  Particles  Adjacent  to 
the  Upper  and  Lower  Wing , Surf aces . 


FIG.  5 CONSTANT  TEMPERATURE  HOT-WIRE  ANEMOMETER 
AND  SQUARE  WAVE  GENERATOR 
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FIG.  9 X-METER  A.C.  SIGNAL  CONDITIONER 


FIG.  10  DETAIL,  X-METER  CIRCUIT  COMPONENTS 


FIG.  11  VARIABLE  ANGLE  X-WIRE  PROBE 
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Fig.  13.  Lineari zed  Hot-Wire 
Response  Characteristic. 
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FIG.  14  TRANSVERSE  VELOCITY  CALIBRATION 
X-METER,  TYPICAL 
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Fig.  15.  Variation  of  Lift  Coefficient  with 
Angle  of  Attack  (NACA  64009) 
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Variation  of  Lift-Drag  Ratio  with  Angle 
of  Attack  (NACA  64009) 


Boundary  layer  thickness  less  than  0.03”  - lower  surface 


O 

> 


0 


4) 

O 

3 

0 


3 


3 

3 

P 

i-i 

rt 

to 

O 

tfl 

o 

0 

I-I 

3 

0 

O 

3 0 

•H 

0 3 

X i — 1 

P 

S 3 

> 

P 

0 

3 

u E 

0 

rt  rt 

G 

4a  0 

O 

3 3 

O 

3 P 

rt 

10  to 

rH 

cx 

E 0 

v> 

O 0 

•H 

3 3 

TJ 

tp  tp 

3 

0 0 

O 

u o 

3 -H 

T) 

3 2 

0 

P P 

(0 

10 

rt 

•H  (A 

pa 

O -n 

O < 


25.  Chordwise  Distribution  of  Boundary  Layer  Thickness 
Inboard  Spanwise  Locations 
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Fig.  28.  Downwasli  Angle  Variation  with  Distance 
from  Upper  Surface,  x/c  = .98 
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Fig.  32.  Crossflow  Velocity  Variation  i 
Inboard  Region,  Lower  Surface 


Fig.  35.  Crossflow  Velocity  Variation  in  the  Inboard 

Region,  Lower  Surface 


Fig.  36.  Crossflow  Velocity  Variation  in  the  Tip  Region 
x/c  = 0.98,  Upper  Surface 


FIG.  37  FLOW  FIELD  NEAR  THE  WING  TIP,  HELIUM  BUBBLE 
TECHNIQUE,  a = 4°,  UQ  = 60  ft/sec 


Fig.  38.  Sketch  of  Streamlines  in  the 
Ouasi-Trefftz  Plane,  x/c  » .99,  a = 4 
Wing  Tip  Region 
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Fig.  42.  Sketch  of  Streamlines  in  the 
Quasi-Tref ftz  Plane,  x/c  = .20,  a = 4° 
Wing  Tip  Region 
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. Sketch  of  Streamli 
reff tz  Plane  , x/c  = . 
Wing  Tip  Region 
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Fig.  44.  Regions  of  "Reverse  Crossflow", 
Wing  Tip  Region,  a = 4°,  U0  = 60  ‘/sec 


Wing  Tip 


View  of  cross-section  in  the  Trefftz  plane 


Fig.  46.  Idealized  Rectangular  "Slab"  in 
Uniform  Flow 
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Fig. 47.  Streamlines  Over  Upper  Surface, 
(0.1  inches  (z/c  = .012)  nominal) 


Streamlines  Over  Lower  Surface 


FIG.  49  FLOW  OVER  THE  INBOARD  REGION,  UPPER  SURFACE, 
HELIUM  BUBBLE  TECHNIQUE,  cs  = 4°, 

UQ  = 60  ffc/sec 
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Fig.  50.  se  of  Contour  Integration 
Adjacent  .o  the  Wing  Surface 
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Fig.  57.  Crossflow  Velocity  Differential  as 
a Function  of  Spanwise  Position 
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Fig.  59.  Comparison  of  Predicted  and  Measured  Values 
of  Crossflow  Differential  Velocity 
at  the  Trailing  Edge 
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Data  in  rectangular  regions  refers  to  percentage 
of  vorticity  contained  in  that  region.  Upper 
number  is  referenced  to  circulation  of  overall 
tip  contour.  Lower  number  (in  parentheses)  is 
referenced  to  maximum  wing  circulation, 


T0  = 5.781  ft2/sec. 


Fig.  62.  Contour  Integrals  in  the  Tip  Region, 
Trefftz  Plane,  x/c  = 0.99 
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Data  in  rectangular  regions  refers  to  percentage 
of  vorticity  contained  in  that  region.  Upper 
number  is  referenced  to  circulation  of  overall 
tip  contour.  Lower  number  (in  parentheses)  is 
referenced  to  maximum  wing  circulation, 

Vn  = 5.781  ft2/sec. 


Fig.  63.  Contour  Integral  in  the  Tip  Region, 
Trefftz  Plane,  x/c  = 0.90. 
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Fig.  64.  Contour  Integrals  in  the  Tip  Region. 
Trefftz  Plane  • 
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Fig.  65.  Contour  Integrals  in  the  Tip  Region, 

TrefCti  Plane 
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Trefftz  Plane 
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Pig.  69.  Ratio  of  Circulation  of  Flow  on  the 
Upper  Surface  to  Corresponding  Region  on 
the  Lower  Surface 
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Fig.  75.  Behavior  of  Vorticity  Field  Near 
a Rounded  Wing  Tip 
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Fig.  79.  Contour  Integrals  and  the  Bound 
Vorticity  Distribution 
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Fig.  82.  Chordwise  (x)  Component  of  Velocity 
Adjacent  to  the  Wing  Surface, 

Inboard  Region,  y/b  = 0.25 
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Fig.  83.  Mainflow  Differential  Velocity  as  a 
Function  of  Chordwise  Location 
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Fig.  84.  Chordwise  Variation  of  Circulation  Function 
as  Obtained  by  "A"  - type  Contours 
(referenced  to  local  spanwise  maximum  value) 
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Fig.  87.  Vorticity  Disi 
Edge  Over  the  Ini 


O-  Contour  Integral  (corrected, 
average  value) 

A-  Cylindrical  Symmetry  Assumed  (1/r) 

Linear  fit  to  contour  data  and 

1/r  distribution 
---  Lifting  Line  Theory 

Lifting  Surface  Theory 

— K.D.  Wood  Approximation 

D--  Rom’s*-52^  Prediction 


) 10  11  12 
a (degrees)' 


Fig.  89.  Variation  of  Lift  Coefficient  with 
Angle  of  Attack  - Theory  and  Experiment 


Fig.  90.  Non-dimensional  Circulation  Function  at  Various 
Downstream  Locations,  a = 4 degrees 


Fig.  91.  Non-dimensional  Circulation  Function  at 
Various  Downstream  Locations,  a = 6 degrees 
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Fig.  93.  Variation  of  Circulation  with  Angle  of 
Attack  for  Various  Downstream  Locations 
(Contour  Integrals) 
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Fig.  96.  Non-dimensional  Circulation  Function  at 
Various  Angles-of-Attack , x/c  = 4.00 
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Fig.  98.  Variation  of  Wake  Shape  in  the  Near  Field, 

a = 6 degrees 
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